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ABSTRACT 

This paper presents a study of the integrated Sachs-Wolfe (ISW) effect from the Planck 2015 temperature and polarization data 
release. This secondary cosmic microwave background (CMB) anisotropy caused by the large-scale time-evolving gravitational 
potential is probed from different perspectives. The CMB is cross-correlated with different large-scale structure (LSS) tracers: 
radio sources from the NVSS catalogue; galaxies from the optical SDSS and the infrared WISE surveys; and the Planck 2015 
convergence lensing map. The joint cross-correlation of the CMB with the tracers yields a detection at 4 a where most of the 
signal-to-noise is due to the Planck lensing and the NVSS radio catalogue. In fact, the ISW effect is detected from the Planck 
data only at 3cr (through the ISW-lensing bispectrum), which is similar to the detection level achieved by combining the 
cross-correlation signal coming from all the galaxy catalogues mentioned above. We study the ability of the ISW effect to place 
constraints on the dark-energy parameters; in particular, we show that Qa is detected at more than 3cr. This cross-correlation 
analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not 
permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, the Planck polarization data are 
used to study the anomalously large ISW signal previously reported through the aperture photometry on stacked GMB features 
at the locations of known superclusters and supervoids, which is in conflict with AGDM expectations. We find that the current 
Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck 
lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we 
have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the 
previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding 
uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the large-scale 
structure traced by the 2MASS Photometric redshift survey (2MPZ) by directly inverting the density field into the gravitational 
potential field. 

Key words. Gosmology: observations - cosmic microwave background - large-scale structure of the Universe - dark engery - 
Galaxies: clusters: general - Methods: data analysis 
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1. Introduction 


This paper, one of a set associated with the 2015 release of 
data from the Planck^ mission, describes the detection and 
characterization of the integrated Sachs-Wolfe (ISW) ef¬ 
fect using external (galaxy-survey catalogues) and internal 
{Planck lensing map) large-scale tracers. The 2015 Planck 
data release offers polarization information on the cosmic 
microwave background (CMB) for angular scales smaller 
than 5°. Whenever possible, this polarization information 
is used to improve our characterization of the ISW signal. 

The ISW effect (Sachs & Wolfe 1967; Rees & Sciama 
1968; Martinez-Gonzalez et al. 1990; Sugiyama 1995) is a 
secondary anisotropy in the CMB, which is caused by grav¬ 
itational interaction of CMB photons with the growing cos¬ 
mic large-scale structure (LSS): 


AT 2 

Tcmb Jo ^ dx 


( 1 ) 


Here, the fractional temperature perturbation 0 is given as 
a line of sight integral over the time-evolving potentials ^ 
in the LSS. The integral is expressed in terms of comoving 
distance y, which is related to the scale factor a according 
to da/dy = a‘^H{a)/c^ with the Hubble function H{a) and 
the speed of light c. The integration is extended to the 
surface of last scattering ycMB — 10 Gpc/h corresponding 
to a redshift of z 1100 in a ACDM cosmology. 

The ISW effect measures the rate of growth of gravita¬ 
tional potentials relative to universes with a critical den¬ 
sity of matter through frequency shifts in the photon dis¬ 
tribution. It is measured by cross-correlating with a tracer 
of the LSS, such as a galaxy catalogue or a reconstructed 
weak gravitational lensing map, in order to distinguish it 
from primary CMB anisotropies; this is because gravita¬ 
tional interaction conserves the Planckian shape of the pho¬ 
ton spectrum. The ISW effect is generated at late times 
when the growth of structure is influenced by a cosmologi¬ 
cal constant, dark energy (Crittenden & Turok 1996), mod¬ 
ified gravity (Hu 2002), or spatial curvature (Kamionkowski 
1996). 

The most direct way of detecting the ISW effect is the 
determination of the cross-correlation or the cross-angular 
power spectrum between the CMB temperature and the 
density of tracer objects such as galaxies. In this way, 
the first detection was reported by Boughn & Crittenden 
(2004) which was subsequently refined by many groups on 
the basis of WMAP data, yielding values for the detec¬ 
tion significance in excess of 4 a (e.g., Fosalba et al. 2003; 
Nolta et al. 2004; Corasaniti et al. 2005; Padmanabhan 
et al. 2005; Vielva et al. 2006; Giannantonio et al. 2006; 
Cabre et al. 2007; Rassat et al. 2007; McEwen et al. 2007; 
Giannantonio et al. 2012). Corresponding constraints on 
cosmological parameters were derived for standard models 
with a cosmological constant and for dark energy models 
(e.g., Pietrobon et al. 2006; McEwen et al. 2007; Vielva 
et al. 2006; Giannantonio 2008; Ho et al. 2008; Xia et al. 


* Corresponding author: P. Vielva vielva@ifca.unican.es 
^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by 
two scientific consortia funded by ESA member states and led 
by Principal Investigators from France and Italy, telescope re¬ 
flectors provided through a collaboration between ESA and a 
scientific consortium led and funded by Denmark, and additional 
contributions from NASA (USA). 


2009), as well as for models with modified gravity (e.g., 
Zhao et al. 2010). A Bayesian ISW detection method, which 
estimates the ISW amplitude conditionally to the observed 
LSS, can be expected to provide 10% better signal-to- 
noise ratio compared to a direct CMB-LSS cross-correlation 
study (Erommert et al. 2008), as used traditionally and in 
this psper beacuse of its lower computational complexity. 

In fact, using the ISW signal alone (but fixing the re¬ 
maining cosmological parameters), the dark energy density 
parameter was estimated to be ~ 0.75 with an error 
of about 20% (e.g., Nolta et al. 2004; Vielva et al. 2006; 
Giannantonio et al. 2006), the dark energy equation of state 
parameter was found to be close to w = —1 (e.g., Vielva 
et al. 2006; Giannantonio et al. 2006; Ho et al. 2008), and 
tests on spatial flatness yielded upper limits of a few per¬ 
cent for {Ik (e.g.. Ho et al. 2008; Li & Xia 2010), thus 
confirming the concordance cosmological model. 

The presence of systematics at large angular scales 
in LSS surveys and their possible impact on ISW stud¬ 
ies was first emphasized in Hernandez-Monteagudo (2010) 
and formally addressed in Giannantonio et al. (2012) and 
Hernandez-Monteagudo et al. (2014). The ISW analysis 
with the Planck data release in 2013 (Planck Collaboration 
XIX 2014) was consistent with WMAP results using the 
NVSS radio catalogue and catalogues of tracer objects de¬ 
rived with optical SDSS data, while lowering the claimed 
detection levels to smaller numbers (from > 4 a down 
to around 2.5 <j). In addition, a non-zero correlation be¬ 
tween the reconstructed CMB-lensing map as an LSS tracer 
and the microwave background was reported for the first 
time, using the non-vanishing bispectrum of the CMB 
anisotropies on the relevant scales. The strength of this 
correlation was measured to be 3 a, and provides fur¬ 
ther evidence for a late-time accelerated expansion of the 
Universe, as theoretically shown by Hu & Okamoto (2002) 
and Okamoto & Hu (2003). 

An alternative method for detecting the ISW effect is 
the stacking of CMB fields at the positions of known su¬ 
perstructures; if the ISW effect is associated with regions 
of large density, it should be possible to reduce the noise 
due to primary, uncorrelated CMB anisotropies by superpo¬ 
sition and to reach a reduction inversely proportional to the 
square root of the number of stacked fields. Detections using 
this method range between 2 a and 4 a, based on WMAP 
data (e.g., Granett et al. 2008a; Papal et al. 2010) and on 
Planck data (Planck Collaboration XIX 2014). 

A third application of the ISW effect is the reconstruc¬ 
tion of a large-scale map of projected gravitational poten¬ 
tials (Barreiro et al. 2008). Using the correlation between 
temperature anisotropies and a map of the tracer den¬ 
sity, it is possible to estimate these secondary temperature 
anisotropies directly. 

The purpose of this paper is the measurement of the 
ISW effect with the full Planck 2015 data set and to es¬ 
tablish the corresponding constraints on cosmological pa¬ 
rameters. In principle, including polarization data allows us 
to reduce the error bars in estimating angular cross-power 
spectra (Erommert & EnBlin 2009), and it provides a sepa¬ 
ration of the temperature anisotropies into those correlated 
and uncorrelated with polarization, through which the sec¬ 
ondary nature of the ISW effect can be better investigated. 
Eurthermore, the reconstruction of the weak lensing poten¬ 
tial is improved, and a better template for cross-correlation 
is provided. However, as mentioned above, the current po- 
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larization information provided in the CMB maps of the 
2015 Planck data release is limited to angular scales smaller 
than 5° (more precisely, only multipoles ^ > 20 are kept, 
with a cosine transition in the range 20 < ^ < 40). This 
limits the amount of information on the ISW effect that 
can be obtained from the polarization data, since this sec¬ 
ondary anisotropy is mostly significant on the largest angu¬ 
lar scales. Therefore, in this paper, polarization is not used 
for the CMB cross-correlation with LSS tracers, although 
it is considered in the analysis of the CMB anisotropies 
stacked on the positions of known superstructures. 

The paper is organized as follows. In Sect. 2 we present 
the data used in this work (both for the CMB and the LSS 
tracers). The cross-correlations of these tracers are inves¬ 
tigated in Sect. 3. In Sect. 4 we present the results of the 
stacking analysis using temperature and polarization data. 
The recovery of the ISW anisotropy map is described in 
Sect. 5. Finally, we discuss our main results and their cos¬ 
mological implications in Sect. 6. 

2. Data sets 

In this section we describe the data sets and the simu¬ 
lations used throughout the paper. In Sect. 2.1 we de¬ 
scribe the CMB related data (temperature and polariza¬ 
tion anisotropies), whereas the LSS data sets are discussed 
in Sect. 2.2, including galaxy, cluster and void catalogues 
from redshift and photometric surveys, and the Planck lens- 
ing map. In Sect. 2.3 we explain the specific simulations 
performed to study the CMB-LSS cross-correlation. 

2.1. CMB data 

There are four major Planck foreground-cleaned CMB tem¬ 
perature and polarization maps, namely, the COMMANDER, 
NILC, SEVEM, and SMICA maps, named after their re¬ 
spectively generating component separation methods (see 
Planck Collaboration IX 2015, for details). All these maps 
are used here, in comparison, in order to test the robust¬ 
ness of our results. Together with the common Q and U 
Stokes parameter polarization maps, the Planck 2015 data 
release also provides F^-mode maps based on the four com¬ 
ponent separation methods. In addition, the SEVEM method 
also provides foreground-cleaned CMB maps at specific fre¬ 
quencies, in temperature at 100, 143, and 217 GHz, and in 
polarization at 70, 100, and 143 GHz. 

The Planck 2015 GMB maps are provided at different 
resolutions (Planck Gollaboration IX 2015). In this paper 
we consider two different resolutions, depending on the ap¬ 
plication. First, maps with a HEALPix (Gorski et al. 2005) 
resolution parameter A^side= 64 (FWHM = 160 arcmin) 
are adopted for studying the GMB-LSS cross-correlation 
(Sect. 3) and for recovering the ISW anisotropies (Sect. 5). 
Second, A'side= 512 (FWHM = 20 arcmin) maps are used 
to study the ISW effect through the stacking of GMB maps 
on the positions of known superstructures (Sect. 4). Each 
resolution has an associated set of masks, one for tempera¬ 
ture (called UT78, /sky = 74% at A^side = 512), another for 
Q and U Stokes parameters (called UPB77, /sky = 76 % at 
A^side = 512), and a final one for the E-mode (/sky = 45 % 
at Aside = 512). The /sky parameter indicates the fraction 
of the sky that is retained after masking. 

In addition, there are 1000 simulations associated with 
each delivered map, which allow us to characterize the in¬ 


strumental properties of Planck GMB maps. In the context 
of this work, these simulations are used for the stacking 
analyses in Sect. 4. The other ISW studies require specific 
coherent simulations between the GMB and the LSS trac¬ 
ers. These simulations are described in Sect. 2.3. 

As mentioned in Sect. I, the polarized GMB maps of 
the 2015 release have been high-pass filtered (see Planck 
Gollaboration VH 2015; Planck Gollaboration IX 2015, for 
details). In particular, all the multipoles with ^ < 20 were 
removed, and a cosine transition between 20 < -^ < 40 was 
imposed. Obviously, this high-pass filtering very much lim¬ 
its the usefulness of the polarization information for the 
ISW analyses. More precisely, the expected 16% increase 
of the ISW detection significance by exploiting polarization 
information in the GMB-LSS cross-correlation (Frommert 
& EnBlin 2009) depends, mainly, on the filtered-out scales 
(up to 80% for i < 20, and 90% for i < 40). In addition, 
the approach to derive the E-correlated {Te-c) and the E- 
uncorrelated {Te-u) maps (see below), are based on an E- 
mode map, with a corresponding mask that, as mentioned 
above, covers significantly less sky (45%) than the tempera¬ 
ture one (74%). Therefore, in practice, there is no real gain 
in the signal-to-noise level. Nevertheless, some of the infor¬ 
mation kept at smaller scales may still be useful for partic¬ 
ular analyses such as the stacking of the GMB anisotropies 
on the positions of known superstructures. Eirst, because 
these structures are within the part of the sky covered by 
the Te-c and Te-u maps and, second, because the mul¬ 
tipole range that mainly contributes to the angular scales 
of the stacked profiles corresponds to smaller scales than 
those for the GMB-LSS cross-correlation. Eor that reason, 
the polarization information is not used in the ISW study 
through the correlation of the GMB and the LSS tracers 
(Sect. 3), but it was considered in the stacking analyses 
(Sect. 4). A final study of the ISW effect using full polar¬ 
ization information is expected to be done with the next 
Planck data release. 

The primary GMB temperature anisotropies act effec¬ 
tively as a noise source for the measurement of secondary 
GMB anisotropies by increasing its cosmic variance. This is 
true for the ISW effect, which does not produce a notable 
E-mode polarization. Hence, polarization data permit us to 
identify the part of the primary temperature anisotropies 
that is correlated with the E-mode polarization, and to 
remove it from the maps. The resulting GMB temperature 
map, partly cleaned form primary anisotropies, provides up 
to a 16 % better signal-to-noise ratio for secondary fluctua¬ 
tions (Erommert & EnBlin 2009). To this end, we separate 
the temperature map in two components: an E-correlated 
(Te-c) and an E-uncorrelated (Te-m) part. Eollowing the 
approach of Erommert & EnBlin (2009), we have produced 
these maps from the delivered GMB inputs described above. 
An estimation of the E-correlated temperature anisotropies 
(Te-c) is given, in terms of its spherical harmonic coeffi- 

dents by 

dm"" = (2) 

where the filter W£ is defined by the TE and the EE angular 
power spectra: 

e^te I eTE 

^ 

CfE ^pEE ^^EE^ W 

with Ci, El, and Ni representing the angular power spectra 
of the GMB, residual foregrounds, and noise, respectively. 
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Fig. 1. Planck CMB temperature and polarization anisotropies as provided by the SEVEN component separation method 
at a resolution of A^side=512. From left to right and from top to bottom, the panels show the maps of temperature,the 
£'-mode, and the F^-correlated (Te-c) and F^-uncorrelated (Te-m) temperature maps. The units are Kelvin. 


Hence, the Te-c map is given by: 

^max i 

Te-c in) = Y^ (n) , (4) 

£=0 m =—£ 

with (n) the spherical harmonic functions. The Te-u 
map is build by subtraction: Te-u (^) = T (n) — Te-c (^)- 
The above procedure is performed by applying an apodized 
version of the corresponding masks. In Fig. 1 we show the 
T, Te-c^ and Te-u maps for SEVEN. In practice, the deter¬ 
mination of the filter W£ is not straightforward; although 
the CMB and noise contributions can be obtained di¬ 
rectly from the Planck best-fit cosmological model (Planck 
Collaboration XIII 2015) and the FFP8 simulations (Planck 
Collaboration XII 2015; Planck Collaboration IX 2015), in¬ 
formation about the residual foregrounds (Fi) present in 
the CMB temperature and polarization is also needed. We 
verified that the expected CMB and noise power spectra ac¬ 
count well for the observed TE and the EE angular power 
spectra at < 200. Although the foreground spectra are not 
fully known, their impact is minor on these scales due to 
the large mask imposed on the F^-mode map and the high- 
pass filtering applied to the polarization data. However, at 
smaller angular scales some foreground residuals exist. 

An alternative way to construct such a filter to reduce 
primary anisotropies is to extract the relevant correlation 
functions directly from the data. In particular, we have 


constructed filters W£ using a smooth fit of the filter con¬ 
structed as the ratio of the TE and the EE angular power 
spectra of the different CMB component separation maps. 
The procedure followed to build the filter distinguishes be¬ 
tween high- and \ow-£ regimes. For small scales {£ > 200), 
we compute the ratio of Cj^ and Cf^ obtained from the 
data using an apodized mask, which is afterwards smoothed 
following the Savitzky-Golay procedure (Savitzky & Golay 
1964). In the low-.^ regime {£ < 200) the filter is constructed 
using the average value obtained from 1000 simulations of 
GMB plus noise, using the same apodized mask. The result¬ 
ing filters (solid lines) are shown in Fig. 2; for comparison, 
the corresponding theoretical filters, computed only from 
the instrumental properties and the Planck fiducial angu¬ 
lar power spectra, are also plotted (dashed lines). 


2.2. LSS tracers 

As mentioned in Sect. 1, tracers of the gravitational poten¬ 
tial of the LSS are required to extract the secondary ISW 
anisotropies from the dominant primary GMB anisotropies. 
These tracers are used to perform the GMB-LSS cross¬ 
correlation, but also for studying the ISW effect through 
the stacking of the GMB anisotropies on the position of 
known superstructures (such as clusters or voids), and for 
producing a map of the ISW anisotropies. 
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We have included three additional galaxy catalogues 
with respect to to the ones used in Planck Collaboration 
XIX (2014), which were the radio NVSS catalogue and 
the optical luminous galaxies (SDSS-CMASS/LOWZ) as 
well as the main photometric galaxy sample (SDSS-MphG) 
catalogues from the Sloan Digital Sky Survey (SDSS). 


These additional catalogues consist of star forming galax¬ 
ies (WISE-GAL), of active galactic nuclei (AGN; WISE- 
AGN), both sets taken from the catalogue of extragalac- 
tic sources detected by the Wide-Field Infrared Survey 
Explorer (WISE, see Wright et al. 2010), and of photomet¬ 
ric redshifts (2MPZ) obtained from the Two Mieron All Sky 



Fig. 2. Eilter used to construct the T^-ccorrelated maps 
for all component separation methods: COMMANDER in red; 
NILC in orange; SEVEM in green; and SMICA in blue. The 
solid lines are obtained directly from the data, whereas the 
dashed ones represent the theoretical shape of the filters, 
only considering the instrumental noise characteristics of 
the data and the fiducial Planek angular power spectra. 

o 



Redshift z 


Fig. 3. Redshift distributions of the different surveys used 
as LSS tracers. To facilitate comparison, the distributions 
of the external tracers have been normalized to unity (and 
multiplied by a factor 10 for the 2MPZ catalogues). For 
completeness, we also include the contribution of the grav¬ 
itational potential to the lensing convergence map, as a 
function of redshift (without any additional normalization). 
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Survey Extended Source Catalogue (2MASS-XSC), WISE 
and SuperCOSMOS data sets. This last catalogue is only 
used to build an estimation of the ISW anisotropies based 
on a reconstruction of the gravitational potential from the 
3D distribution of the galaxies (see Appendix. A). The rea¬ 
son is that the expected CMB-LSS cross-correlation sig¬ 
nal is very low to be used in this cross-correlation study 
but, however, the galaxy redshift estimation error is suf¬ 
ficiently low to attempt the gravitational potential recon¬ 
struction. Einally, we also cross-correlate the Planck lensing 
map as a LSS tracer with the CMB. In particular, we use 
the lensing convergence map (Kappa) obtained in Planck 
Collaboration XV (2015). 

The redshift distributions of these catalogues are shown 
in Eig. 3. We note that lensing, NVSS, and WISE-AGN 
offer the widest redshift coverage. Some basic properties 
of the galaxy catalogues used (NVSS, WISE-AGN, WISE- 
GAL, SDSS-GMASS/LOWZ, SDSS-MphG, and 2MPZ) are 
summarized in Table 1. 

Eor a better visualization, Wiener-filtered versions of 
the all-sky density projection of the external catalogues, 
as well as the Planck Kappa map, are shown in Eig. 4. 
These are constructed from the theoretical power spectra 
obtained as described in Sect. 2.3. In Eig. 5, we show the 
angular auto- and the cross-power angular spectra for all 
the LSS tracers: dashed lines and points correspond to the 
theoretical model and the data measurements, respectively 
(red for auto-spectra, and blue for cross-spectra); and grey 
areas represent the 1 a sampling uncertainties due to cos¬ 
mic variance. All of these spectra have been corrected for 
the mask coupling following the MASTER approach (Hivon 
et al. 2002). Notice that the Planck lensing convergence 
map only contains information for multipoles £ > S (see 
Planck Gollaboration XV 2015, for details). The two maps 
based on the WISE catalogues (WISE-AGN and WISE- 
GAL) exhibit some extra signal at the largest scales. We 
identify this with some systematic effect present in these 
catalogues and, therefore, as a baseline, we only consider 
multipoles ^ > 9 for these two surveys. This cut implies 
only a minor loss of the ISW signal, while permitting a more 
robust determination of it. The rest of the auto-spectra are 
in reasonably good agreement with the theoretical predic¬ 
tions. Notice that any mismatch on the auto-spectra could 
suffer, not only from systematic effects, but from an inac¬ 
curate description of the statistical properties of the cat¬ 
alogues. In this sense, cross-spectra are, in principle, less 
affected by systematics (at least, among catalogues from 
different experiments), and, therefore, are more useful for 
identifying possible problems in the adequacy of the galaxy 
redshift distribution and galaxy biases. We emphasise that, 
in this sense, the Planck Kappa map could in principle be 
a more robust LSS probe, since it does not suffer from 
these kinds of uncertainty and, therefore, its correlation 
with the rest of the surveys is very useful for highlighting 
potential issues related to the catalogue characterization. In 
this sense, from Eig. 4, it seems that the measured cross¬ 
correlation of the lensing potential with the galaxy cata¬ 
logues is very good, indicating that, within the current un¬ 
certainties, the description of the surveys is accurate. The 
three maps (Kappa, WISE-AGN, and WISE-GAL) shown 
in Eig. 4 do not include the cut multipoles. 

Besides the galaxy surveys described above, we also use 
superstructure catalogues to study the ISW effect through 
the stacking of the GMB anisotropies on the positions of 


clusters and voids. We concentrate on the supercluster and 
void catalogue of Granett et al. (2008b), obtained from 
SDSS (GR0808), since, as shown in Planck Gollaboration 
XIX (2014), its reported strong signal would be a challenge 
for the standard AGDM cosmology if it is solely caused by 
the ISW effect. 

Below we provide a description of all these LSS trac¬ 
ers. Eor those catalogues already used in our previous pub¬ 
lication (NVSS, SDSS-GMASS/LOWZ, SDSS-MphG, and 
GR0808) only a summary is provided; more detailed de¬ 
scription can be found in Planck Collaboration XIX (2014). 


2.2.1. The NVSS radio-galaxies catalogue 


The luminous AGN are very powerful radio sources, which 
can be seen also at high redshifts. These sources are able to 
trace the cosmic density field for both the redshift evolu¬ 
tion and the spatial distribution. Therefore they can probe 
the spatial distribution of large-scale potential wells that 
contribute to the ISW effect during the dark energy era. 

The sources we use in this paper are the same sam¬ 
ples we used in Planck Collaboration XIX (2014), i.e., 
the NRAO Vary-Large-Array (VLA) Sky Survey (Condon 
et al. 1998). This NVSS survey was conducted by using 
VLA at 1.4 GHz, and covers up to an equatorial latitude of 
6 e = —40°, with an average noise level of 0.45 mJy beam“^. 
There are roughly 1.4 x 10® sources above a flux thresh¬ 
old of 2.5 mJy. Eigure 4 (top-left panel) shows the all-sky 
density projection for the NVSS galaxies, where the grey 
area indicates regions not observed or disregarded by the 
surveys. Eigure 5 includes a subplot to show the angular 
power spectra (blue points) of the NVSS survey. 

Eor the galaxy bias, we use the Gaussian bias evolution 
model of Xia et al. (2011), i.e., the bias of the survey is 
given by a mass-weighted average. 


b{z) = 


b{M,z) n{M,z) 


( 5 ) 


where n(M, z) is the halo mass function for which we 
adopt the Sheth-Tormen (Sheth & Tormen 1999) form and 
5(M, z) is the bias of halos with comoving mass M. This 
bias (as a function of redshift) can be approximated as a 
second-order polynomial, as given in Table 1. In addition, 
the redshift distribution is parametrized by: 




(6) 


where zq = 0.33 and a = 0.37, and no is a constant to 
normalize the distribution to unity. The function is given by 
the red line in Eig. 3. We refer the interested readers to our 
previous paper Planck Collaboration XIX (2014) for more 
details of the possible systematic effects for the samples. 


2.2.2. The Sloan Digital Sky Survey catalogues 

We use two subsamples of SDSS: the SDSS luminous galaxy 
samples (SDSS-GMASS/LOWZ); and the main photomet¬ 
ric SDSS galaxy sample (SDSS-MphG). The redshift dis¬ 
tributions of the two samples are shown in blue and green 
lines in Eig. 3. The sky coverage for the two subsamples 
are shown in the subplots of Eig. 4, and the angular power 
spectra in the subplots of Eig. 5. 
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Fig. 4. Density contrast maps obtained from the galaxy catalogues at A^side = 64. From left to right and from top to 
bottom: NVSS; WISE-AGN; WISE-GAL; SDSS-CMASS/LOWZ; SDSS-MphG; and 2MPZ. The Planck lensing conver¬ 
gence map (Kappa) is given in the fourth row. Eor visualization purposes, all these maps are Wiener-filtered versions of 
the original data. Maps are in dimensionless units here. 
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Table 1. Main characteristics of the galaxy catalogues used as tracers of the gravitational potential. From left to right, the 
columns indicate the number of galaxies per steradian, the fraction of the sky covered by each survey, the galaxy redshift 
distribution (dn/dz\ the galaxy bias, and the mean redshift. Whereas for the NVSS and the SDSS-MphG catalogues there 
are analytical expressions of the galaxy redshift distribution, for the other tracers there are only numerical estimations. 


Galaxy catalogue 

n 


/sky 

dnjdz 

bias 

Z 

NVSS 

1.584 

X 


0.73 

oc(z/0.33)°“e-o®"("/°®®) 

0.90 (1 +0.54(1 + z)^) 

1.22 

WISE-AGN 

1.552 

X 


0.45 

numerical 

0.48 (1 +0.54(1 + z)/ 

1.03 

WISE-GAL 

1.187 

X 

10 ^ 

0.45 

numerical 

0.79 (l + «) 

0.38 

SDSS-CMASS/LOWZ 

5.558 

X 


0.22 

numerical 

2.03 

0.45 

SDSS-MphG 

9.680 

X 

10 ® 

0.22 


1.20 

0.32 

2MPZ 

8.328 

X 

10"^ 

0.66 

numerical 

1.35 

0.09 


SDSS-CMASS/LOWZ 

We use the photometric luminous galaxy (LG) cata¬ 
logue from the Baryonic Oscillation Spectroscopic Survey 
(BOSS) of the SDSS III. The data used consist of two sub¬ 
samples: GMASS; and LOWZ. In this paper we will use 
a combination of them, i.e., SDSS-GMASS/LOWZ in our 
data analysis. 

The GMASS sample has roughly constant stellar mass, 
and is mostly contained in the redshift range 2 : = 0.4-0.7, 
with a galaxy number density close to 110deg~^. With the 
colour selection criteria, it is a catalogue of about one mil¬ 
lion sources, in an area of 10500 deg^. Photometric redshifts 
of this sample are calibrated using a selection of about 
100,000 BOSS spectra as a training sample for the pho¬ 
tometric catalogue. 

The photometric LOWZ sample selectes luminous, 
highly biased, mostly red galaxies, placed at an average 
redshift of z ^ 0.3 and below the redshifts of the GMASS 
sample {z < 0.4). With a total number of sources of roughly 
600000, the number density of galaxies in the southern part 
of the footprint is higher than in the northern one (by more 
than 3%). Both SDSS-GMASS and SDSS-LOWZ samples 
are further corrected for any scaling introduced by possible 
systematics like stars, mask value, seeing, sky emission, air 
mass and dust extinction, since the high star density tends 
to “blind” galaxy detection algorithms. The algorithm fol¬ 
lowed to correct for systematics is described in Hernandez- 
Monteagudo et al. (2014). 

The LOWZ and GMASS galaxy samples come from two 
different colour and magnitude selections on SDSS pho¬ 
tometry, and they effectively probe two different redshift 
ranges. If studied individually, the expected sensitivity of 
GMASS is larger than for LOWZ, resulting in the total 
signal-to-noise obtained when adding the contribution of 
each survey separately being very similar (to 10-15%) to 
the result from a combination of the two surveys into one 
single galaxy sample. Most of this difference comes from 
assuming an effectively constant bias for the joint survey. 
Therefore, the improvement is quite low, particularly tak¬ 
ing into account that this is a survey that, at the end of the 
day, provides very little of the signal-to-noise in the total 
4cr detection (see Sect. 3). On the other hand, the effective 
combination simplifies the overall analysis. Regarding the 
mask, we also remark that, since it is determined by sys¬ 
tematic effects such as the star density, airmass, or Galactic 
extinction, and is built independently of the colour and 
magnitude cuts applied on the galaxy sample from SDSS, 
we ended up having the same mask for both the LOWZ 


and GMASS samples (even if one could argue for more or 
less conservative masks). 

SDSS-MphG 

These are the photometrically selected galaxies from the 
SDSS-DR8 catalogue, which covers a total sky area of 
14 555deg^ (Aihara et al. 2011). The total number of ob¬ 
jects labelled as galaxies in this data release is 208 million. 
However, for correcting extinction and restricting redshift 
ranges, our final sample consists of about 42 million, with 
redshifts distributed around a median value of 0.35. We 
use the analytical function 

dz Y (^2^) 

with parameters m = 1.5, [3 = 2.3, and zq = 0.34, for 
the number density distribution, and the constant galaxy 
bias b = 1.2 by fitting the AGDM prediction to the ob¬ 
served auto-correlation function of the galaxies. As for 
the LOWZ/GMASS LRG samples, this galaxy sample was 
also corrected for systematics following the approach of 
Hernandez-Monteagudo et al. (2014). 


2.2.3. The Wide-Field Infrared Survey Explorer extragalactic 
catalogues 

We next describe the use of the extragalactic sources de¬ 
tected by the textitWide-Field Infrared Survey Explorer 
(WISE, see Wright et al. 2010) in our ISW studies. WISE 
scanned the full sky at 3.4, 4.6, 12 and 22 pm (which con¬ 
stitute bands W1 to W4). These observations give a deeper 
view of the infrared sky than previous surveys like 2MASS 
or IRAS, and provide an extensive extragalactic catalogue 
of more than 500 million sources (see Wright et al. 2010). 
The four W1 to W4 bands are sensitive to either UV ra¬ 
diation reproduced by dust grains in star forming galaxies 
or to infrared emission from stars, either in our Galaxy or 
in extragalactic sources. The W1 band turns out to be the 
deepest, sampling the deep Universe by detecting massive 
galaxies up to z 1 and with a median redshift of 0.3 (Yan 
et al. 2013) 

In the context of ISW studies, our approach is very sim¬ 
ilar to that of Eerraro et al. (2014). We focus our efforts 
on two different sets of extragalactic sources: star-forming 
galaxies; and AGN. This requires a careful separation of the 
stars in the catalogue. Since the sky scanning of the WISE 
satellite is not homogeneous, a magnitude cut of Wl< 16.6 
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Fig. 5. Angular power spectra from the maps in Fig. 4 used to study the ISW effect through the CMB-LSS cross¬ 
correlation. From top to bottom: NVSS; SDSS-CMASS/LOWZ; SDSS-MphG; WISE-GAL; WISE-AGN; 2MPZ; and 
Kappa. The observed spectra are the points (red for auto-spectra and blue for the cross-spectra), while the theoretical 
models are represented by the dashed lines (the grey areas correspond to the sampling variance). 


is imposed on all sources at high galactic latitude, since 
for this cut Eerraro et al. (2014) found a uniform sample. 
Given that stray light from the Moon may cause faint de¬ 
tections and other spurious effects in the data, we discard 
all sources with a flat moon_lev > 4, while also dropping 
all sources suspected of being artefacts (cc_f lags?^ 0). 

Eollowing the colour cuts given in Yan et al. (2013) and 
Eerraro et al. (2014), we impose the cut Wl-W2> 0 to 
isolate the galaxies from stars, while the stricter conditions 
Wl—W2> 0.85 and W2<15.0 should separate the AGN 
from star forming galaxies, although the former constitute a 
very small fraction of the latter. With these cuts, we obtain 
about 140 million galaxies and 1.4 million AGN in the entire 
sky. 

The presence of systematics causes clear excess in the 
auto-power spectra of these two WISE-based surveys on 
the largest scales. This poses a problem, since, in order to 
predict the level of cross-correlation with the GMB maps 
generated by the ISW component, we need first to char¬ 
acterise the bias of each tracer. Eor this purpose we follow 
exactly the same approach as in Eerraro et al. (2014): we use 


the cross-correlation of these two galaxy surveys with lens- 
ing convergence maps from Planck in order to estimate the 
bias. This requires adopting some models for the redshift 
distribution of WISE galaxies and AGN, and in particular 
we use those given in Yan et al. (2013), which were obtained 
after cross-matching WISE and SDSS data. We also adopt 
some redshift dependences for the bias that are identical 
to those used in Eerraro et al. (2014). Eor WISE-GAL our 
fiducial model has this very simple redshift dependence: 

^WISE-GAL ^ j,WISE-GAL ^ (g) 

Eor the WISE-AGN the suggested redshift dependence is 
quadratic: 

^WISE-AGN yj ^ ^WISE-AGN [q 53 ^ g 289(1 + ^)^]. (9) 

After cross-correlating with Planck lensing convergence 
maps in the multipole range £ = 10—400, we find the follow¬ 
ing values for the fiducial parameters: ^wise-gal ^ q yg. 
and ^ 0.90. 
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2.2.4. The 2MASS Photometric Redshift catalogue 

The 2MASS Photometric Redshift catalogue (hereafter 
2MPZ Bilicki et al. 2014) is a combination of the 2MASS 
XSC (Jarrett et al. 2000), WISE, and SuperCOSMOS 
(Hambly et al. 2001) surveys, yielding an all sky extragalac- 
tic source catalogue with a typical uncertainty in redshift of 
az = 0.016. This is achieved by employing an artificial neu¬ 
ral network approach (the ANNz algorithm, see Collister & 
Lahav 2004) on the above quoted surveys, and after train¬ 
ing it with the 2MRS (Huchra et al. 2012), SDSS^, 6dFGS 
(Jones et al. 2009), 2dFGRS^, and ZGAT (Huchra et al. 
1995) spectroscopic surveys. 

The resulting catalogue contains almost one million 
sources with a median redshift of 0.09. Out of those sources, 
more than three hundred thousand contain spectroscopic 
redshifts. In Fig. 3 we display the histogram of the pho¬ 
tometric redshift distribution of these sources, with a high 
redshift tail extending up to z 0.3. When using MASTER 
(Hivon et al. 2002) to compute its angular power spectrum, 
we find that its effective linear bias for large scales {£ < 70) 
is close to unity (6 1.35). However, at small scales there 

is clear evidence for non linear power. 

2.2.5. Planck lensing map 

The clustered matter forming the cosmic web modifies the 
2D projected distribution of the GMB anisotropies on the 
sky, via the weak gravitational lensing effect. This distor¬ 
tion breaks the isotropy of the intrinsic GMB fluctuations, 
introducing correlations among multipoles. Optimal inver¬ 
sion methods (Hu & Okamoto 2002; Okamoto & Hu 2003) 
allows us to recover the projected density field (0), which 
is proportional to the gravitational field (4>). 

As part of its official products release, Planck provides 
a map of the estimated lensing field (Planck Collaboration 
XVH 2014; Planck Collaboration XV 2015), which can be 
used to probe the ISW effect, through its cross-correlation 
with the GMB map, in the same manner as is done with ex¬ 
ternal galaxy catalogues. The lensing signal measured from 
this map is detected at about 40 a, using the full-mission 
temperature and polarization data to construct the lensing 
map estimation. 

In fact, what Planck releases is the lensing convergence 
map (/^), which has a whiter angular power spectrum than 
the raw lensing potential (0): Kirn = + l)/2- This 

lensing map was obtained from the SMICA GMB solution, 
and it is shown (Wiener-filtered) in Fig. 4, bottom row. It 
covers 67% of the sky, and has a multipole range S < £ < 
2048. The convergence map traces the matter distribution 
through a wide redshift range (see Fig. 3). 

2.2.6. The Sloan Digital Sky Survey superstructures 

We use here the catalogue of superstructures^ from Granett 
et al. (2008a), also used in Planck Collaboration XIX 
(2014). This sample consists of 50 superclusters and 50 
supervoids identified from the Luminous Red Galaxies 
(LRGs) in the SDSS (sixth data release, DR6, Adelman- 
McCarthy et al. 2008), which covers an area of 7500 deg^ on 

^ http://www.sdss.org 
^ http://www2.aao.gov.au/~TDFgg/ 

^ http://ifa.hawaii.edu/cosmowave/supervoids/ 


the sky. These authors used publicly available algorithms, 
based on the Voronoi tessellation, to find 2836 super¬ 
clusters (using VOBOZ, VOronoi BOund Zones, Neyrinck 
et al. 2005) and 631 supervoids (using ZOBOV, ZOnes 
Bordering On Voidness, Neyrinck 2008) above a 2 a sig¬ 
nificance level (defined as the probability of obtaining, in 
a uniform Poissonian point sample, the same density con¬ 
trasts as those of clusters and voids). 

The 50 superclusters and 50 supervoids published in 
the Granett et al. (2008a) catalogue correspond to density 
contrasts of about 3 a and 3.3 a, respectively. They span a 
redshift range of 0.4 < ^ < 0.75, with a median of around 
0.5, and inhabit a volume of about 1.6 Gpc^. These super¬ 
structures can potentially produce measurable ISW signals, 
as suggested in Granett et al. (2008a, b). For each structure, 
the catalogue provides: the position on the sky of its cen¬ 
tre; the mean and maximum angular distance between the 
galaxies in the structure and its centre; the physical vol¬ 
ume; and three different measures of the density contrast 
(calculated from all its Voronoi cells, from only its over- or 
under-dense cells, and from only its most over- or under- 
dense cell). 

2.3. Simulations 

We have performed a set of 11000 correlated simulations 
of the GMB and different LSS tracers, which are used to 
study the GMB-LSS cross-correlation up to ^ ~ 190. In fact, 
the maps are produced directly at A^side= 64, since we have 
checked that these maps already capture all the information 
required for our analyses. From the GMB side, each sim¬ 
ulation consists of two independent signals: an ISW map; 
and the rest of GMB anisotropies in a T map. Notice that, 
since no polarization information at the largest scales is pro¬ 
vided in this release, this signal is not generated. At these 
large scales, instrumental noise is negligible and, there¬ 
fore, we have not included it in the simulations. From the 
galaxy-surveys side, we simulate galaxy density maps for 
NVSS, SDSS-GMASS/LOWZ, SDSS-MphG, WISE-AGN, 
and WISE-GAL. Although it is not used for studying the 
GMB-LSS cross-correlation, we also generate a 2MPZ map 
that helps to assess the quality of the ISW recovery from 
this photometric catalogue (see Sect. 5). Shot noise is added 
to the simulated galaxy density maps, accordingly to the 
corresponding mean number of galaxies per pixel given in 
Table 1. Finally, we also produce a coherent lensing conver¬ 
gence (Kappa) map. Noise is also added to this map, fol¬ 
lowing the uncertainty level estimated in the Planck lensing 
paper (Planck Gollaboration XV 2015). 

The maps are simulated by assuming that the pertur¬ 
bations are Gaussian. This is a good approximation for 
the Kappa map, as well as for all the distributed galaxy 
density maps. In particular, for the galaxy catalogues, al¬ 
though they follow a Poisson distribution, the mean number 
of galaxies per pixel is large enough (about 40, for the worst 
case, 2MPZ). Therefore, all the required information to per¬ 
form the coherent simulations is given by all the angular 
auto- and cross- power spectra (see, for instance, Barreiro 
et al. 2008, for details). In particular, given two surveys 
a and 6, the theoretical cross-power spectra between the 
surveys read as: 

Cf = 47r I ^ A\k)ink)lKk) , (10) 
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where A‘^{k) is the matter power spectrum per logarithmic 
interval, and I^{k) is a transfer function represented by the 
redshift integral: 

nOO 

It{k)= / dzWa{z,k)j^{kr{z)) . ( 11 ) 

Jo 

Here r{z) is the comoving distance as a function of the 
redshift, and are the spherical Bessel functions, which 
project the window function Wa{z, k) into each multipole I 
of the power spectrum. In the case of a galaxy survey, the 
window function is independent of k and is given by 

dn 

Wa(z) = ba{z)D+(z)^ . ( 12 ) 

This depends on the galaxy redshift distribution and the 
bias function ba{z) of each survey. The growth factor D_^{z) 
in this expression takes into account the linear evolution of 
the matter perturbations. 

For leasing, the efficiency window function H4:(x), 
which relates the density perturbations 6 to the weak leas¬ 
ing convergence Kappa in a line of sight integration, 

K = JdxW^{x)5, (13) 

is given by 

= — — —(XCMB-x), (14) 

2Xff a XCMB 

with the comoving distance Xcmb to the surface of last 
scattering, which is approximately 15 Gpc in the ACDM 
cosmology considered here and xh = c/Hq is again the 
Hubble distance. In order to compare the leasing efficiency 
function W^{x) to the redshift distributions of other LSS 
surveys, we convert it by analogy into a dimensionless func¬ 
tion by multiplying it by the inverse Hubble parameter: 

W,iz)dz = W,{x)dx ^ W,{z) = W^.(X)^ = 

(l4 

For the ISW effect, the window function involves the 
evolution of the potential with redshift: 

Wisw{z,k) =■^[{1 +z)D+{z)] . (16) 

This depends on k due to the Poisson equation relating 
the matter and the potential. If the Universe is matter- 
dominated, then the function (1 P z)D_^{z) is constant and 
the ISW effect vanishes. 

All the angular power spectra used in the present pa¬ 
per have been calculated using a modified version of the 
CAMB^ code. The fiducial ACDM cosmological model as¬ 
sumed is: = 0.0222, = 0.119, Vtyh? = 0, 

VLk = 0, ris = 0.9615, As = 2.1740 x 10“^ r = 0.077, 
and h = Hq/100 kms“^Mpc“^, fully compatible with the 
Planck fiducial model (Planck Collaboration XHI 2015). 


3. CMB correlation with tracers of the 
gravitational potential 

The CMB cross-correlation with LSS tracers of the mat¬ 
ter distribution is the most classical approach to study 
the IWS effect; it offers the possibility of extracting these 
secondary anisotropies, otherwise swamped by the primor¬ 
dial CMB anisotropies. The seminal work by Crittenden 
& Turok (1996) proposed to use a galaxy catalogue as the 
LSS tracer, and the first positive detection following this 
approach was done by Boughn & Crittenden (2004) using 
WMAP data and radio and X-ray tracers. As discussed in 
Sect. 1, several other works quickly came after this, confirm¬ 
ing the detection of the ISW effect with additional galaxy 
tracers. 

In our previous study (Planck Collaboration XIX 2014), 
we performed, for the first time, the detection of the ISW ef¬ 
fect using only CMB data, by cross-correlating the Planck 
CMB map with the Planck leasing potential map, which 
is naturally used as a tracer of the matter distribution. 
This cross-correlation is nothing but an estimator of the 
ISW-lensing bispectrum induced on the Caussian CMB 
anisotropies by the deflection caused by the leasing ef¬ 
fect (see, e.g., Lewis et al. 2011). 

This cross-correlation has been studied using different 
tools: the cross-correlation function (CCF); the covariance 
of wavelet coefficients (Wcov); and the angular cross-power 
spectrum (ACPS). These three estimators were first com¬ 
pared in Vielva et al. (2006), and were also applied to study 
the cross-correlation of the Planck CMB data with several 
surveys in Planck Collaboration XIX (2014). 

The ACPS is the most natural tool for studying the 
cross-correlation, since the entire cross-correlation signal is 
fully included in it. In theory, the CCF and the Wcov esti¬ 
mators are also optimal, as long as they are evaluated at a 
sufficient number of angles/scales. In fact, for a given case, 
the CFF and the Wcov, using a relatively small number of 
evaluations, can achieve a significant fraction of the total 
signal-to-noise. The clear advantage of ACPS is that, under 
certain conditions, it provides statistics with uncorrelated 
elements in cases that full-sky maps can be used. Even for 
incomplete sky signals, an angular pseudo-spectrum, ob¬ 
tained through a MASTER approach (e.g., Hivon et al. 2002; 
Hinshaw et al. 2003), provides a very good approximation 
to our problem. 

We showed in Planck Collaboration XIX (2014) that 
the three estimators (CCF, Wcov, and ACPS) render a 
similar detection level, although ACPS obtained tighter 
limits than the other two approaches (especially as com¬ 
pared to the CCF). For that reason, in this release we 
only use the ACPS estimator. The measured ACPS be¬ 
tween SEVEM and the LSS tracers considered in this analysis 
(NVSS, WISE-ACN, WISE-CAL, SDSS-CMASS/LOWZ, 
SDSS-MphC, and Kappa) are shown in Eig. 6. 

3.1. Methodology 

We aim to study the ISW cross-correlation by estimating 
the best-fit amplitude of the ACPS for a given fiducial 
model (the one mentioned in Sect. 2.3). This approach al¬ 
lows us to check the compatibility of the data with the ISW 
effect, and provides an estimate of the signal-to-noise ratio 

^ http://camb.info 
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Fig. 6. Measured ISW-LSS cross-spectra (CAPS). From left to right, snd top to bottom, the panels show the cross¬ 
correlation of the four CMB maps with NVSS, WISE-AGN, WISE-GAL, SDSS-GMASS/LOWZ, SDSS-MphG, and 
Kappa. Grey areas represent ±la uncertainities derived from simulations. Spectra derived from the Planck GMB maps 
are virtually the same. 


of the measured signal. This method is complementary to 
an alternative approach, in which the signal is compared 
to the null hypothesis of no correlation between the GMB 
and the LSS. Using a Bayesian hypothesis test, we showed 
in Planck Gollaboration XIX (2014) clear evidence in favour 
of the alternative hypothesis as compared to the null one. 

Let us denote the expected AGPS of two maps {x and y) 
by where £ represents a given multipole, and we assume 
that the two signals are given in terms of a fluctuation field 
(i.e., with zero mean and dimensionless units). 

In our particular case, x can be seen as the GMB signal, 
and y represents for one or more surveys. In other words, 
we can pursue the estimation of the ISW amplitude by a 
single correlation of the GMB with a given survey, or with 
several surveys jointly. In this latter case, is a vector of 
^max components, where the first ^maxi components corre¬ 
spond to the GMB cross-correlation with the first survey, 
the next ^max 2 components correspond to the correlation 
with the second survey, and so on. Obviously, when x = y, 
represents an auto-correlation. 

The full description of and its covariance in 

terms of the theoretical model and the specific sky cov¬ 
erage, is given in Planck Gollaboration XIX (2014). Since 
the AGPS is a very fast estimator, in particular for the 
Aside = 64 resolution maps, it is also possible to determine 
these quantities from the coherent simulations described in 
Sect. 2.3: we use 10 000 out of the 11000 performed simu¬ 
lations to estimate both the expected signal (C^^) and its 
covariance (C^xy). This is the approach followed in the past 
Planck release, and it is also the one adopted in this work. 

Denoting the observed cross-correlation by a simple 

can be formed to estimate the amplitude A, such that 


A X is the best-fit solution to q 

1 T 


xy. 


X^{A)= er 




xy 


— 1 


£xy 


A X 


cxy 


(17) 


where C^xy is the covariance matrix (of dimension -^max x 
'^max) of the expected cross-correlation i.e., C^xy_ = 

straightforward to show that the best-fit 
amplitude A, its error, and the significance are given by 


A = 

(JA = 


1 T 


(—1 cxy \\dXy]T .— 1 ^xy 
- 1/2 

L (^-^y 

£x 


icxyA" (—1 cxy 


(18) 
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3.2. Cross-correlation results 

The fundamental GMB-LSS cross-correlation results are 
summarized in Table 2, where we report the estimated 
ISW amplitude (A), its error and the detection level 

Aj GAi derived for the surveys described in 2.2, applying 
Eqs. (18). Results obtained from the four Planck GMB 
maps (COMMANDER, NILC, SEVEM, and SMICA) are given, 
showing perfect agreement among them, indicating a ro¬ 
bust recovery of the largest GMB anisotropies. 

In addition to estimating the ISW amplitude by fit¬ 
ting individual surveys, we also consider several combina¬ 
tions: NVSS and Kappa; the two SDSS surveys (SDSS- 
GMASS/LOWZ and SDSS-MphG); the two WISE cat¬ 
alogues (WISE-AGN and WISE-GAL); the five external 
tracers (NVSS, and WISE and SDSS surveys); and the six 
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Fig. 7. Dependence of the estimated Aisw amplitude for the different surveys as a function of the ^min considered in 
the amplitude estimation. 

Table 2. ISW amplitudes A, errors a a, and significance levels S/N = A/ a a of the CMB-LSS cross-correlation (survey- 
by-survey and for different combinations). These values are reported for the four Planck CMB maps: COMMANDER; NILC; 
SEVEM; and SMICA. The last column gives the expected S/N within the fiducial ACDM model. 


LSS data 

COMMANDER 


NILC 


SEVEM 


SMICA 


Expected 


A±aA 

S/N 

A±aA 

S/N 

A±aA 

S/N 

A±aA 

S/N 

S/N 

NVSS 

0.95 ±0.36 

2.61 

0.94 ±0.36 

2.59 

0.95 ±0.36 

2.62 

0.95 ±0.36 

2.61 

2.78 

WISE-AGN (£n.in > 9) 

0.95 ±0.60 

1.58 

0.96 ±0.60 

1.59 

0.95 ±0.60 

1.58 

1.00 ±0.60 

1.66 

1.67 

WISE-GAL > 9) 

0.73 ±0.53 

1.37 

0.72 ±0.53 

1.35 

0.74 ±0.53 

1.38 

0.77 ±0.53 

1.44 

1.89 

SDSS-CMASS/LOWZ 

1.37 ±0.56 

2.42 

1.36 ±0.56 

2.40 

1.37 ±0.56 

2.43 

1.37 ±0.56 

2.44 

1.79 

SDSS-MphG 

1.60 ±0.68 

2.34 

1.59 ±0.68 

2.34 

1.61 ±0.68 

2.36 

1.62 ±0.68 

2.38 

1.47 

Kappa (^min > 8) 

1.04 ±0.33 

3.15 

1.04 ±0.33 

3.16 

1.05 ±0.33 

3.17 

1.06 ±0.33 

3.20 

3.03 

NVSS and Kappa 

1.04 ±0.28 

3.79 

1.04 ±0.28 

3.78 

1.05 ±0.28 

3.81 

1.05 ±0.28 

3.81 

3.57 

WISE 

0.84 ±0.45 

1.88 

0.84 ±0.45 

1.88 

0.84 ±0.45 

1.88 

0.88 ±0.45 

1.97 

2.22 

SDSS 

1.49 ±0.55 

2.73 

1.48 ±0.55 

2.70 

1.50 ±0.55 

2.74 

1.50 ±0.55 

2.74 

1.82 

NVSS and WISE and SDSS 

0.89 ±0.31 

2.87 

0.89 ±0.31 

2.87 

0.89 ±0.31 

2.87 

0.90 ±0.31 

2.90 

3.22 

All 

1.00 ±0.25 

4.00 

0.99 ±0.25 

3.96 

1.00 ±0.25 

4.00 

1.00 ±0.25 

4.00 

4.00 


surveys together. As expected, the lowest error is achieved 
by combining all the surveys, taking into account all their 
mutual correlations. For the fiducial ACDM model a total 
signal-to-noise of 4 a is predicted, and, that is the actual 
value estimated from the data. 

The highest contribution comes from the Planck con¬ 
vergence lensing map (Kappa), which provides a detection 
level of 3.2 cr, followed by NVSS that allows us to detect 
the ISW effect at 2.6 a. In fact, the combination of these 


two LSS tracers almost provides the full detection achieved 
with the six surveys, 3.8 a. 

The ISW effect characterized from the SDSS catalogues 
has a signal-to-noise level of around 2.4 a for each sur¬ 
vey, and 2.7 cr when they are considered jointly. The WISE 
surveys provide the lowest signal-to-noise: 1.6 cr for WISE- 
AGN; 1.4 cr for WISE-GAL; and 1.9 cr for the combination 
of both. The signal-to-noise achieved by the combination 
of the five external tracers is 2.9 cr. All these detection lev- 
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els refer to SEVEM, although the levels achieved from the 
analysis of the other Planck CMB maps are virtually the 
same. 

All the estimated amplitudes are compatible with unity, 
within the corresponding Icr level. In fact, the value of the 
ISW amplitude is quite stable, independent of the lowest 
multipole (^min) considered in the amplitude estimation. 
This is graphically represented in Fig. 7, where the best-fit 
amplitude A (solid-blue lines) and the la error (grey areas) 
are shown. It is remarkable that the estimated amplitude 
is very constant as a function of ^min, and compatible with 
unity. The case in which there is a mild incompatibility is 
for WISE-AGN, where, for i^nin ^ 18, the departure from 
unity is in tension at about la. This could be the indication 
of some systematics or contamination, still present in this 
catalogue. Looking at Fig. 5, we see that the WISE-AGN 
auto-spectra presents some deviations with respect to the 
fiducial model, not only at the lowest multipoles, which, as 
already mentioned, are removed from the analyses, but, in 
general, over the whole £ range, showing some extra power. 
Interestingly, there is no systematic discrepancies between 
the cross-spectra of WISE-AGN and the rest of the cat¬ 
alogues. The same is observed in the cross-correlation of 
GMB with WISE-AGN in Eig. 6. This could indicate that 
the WISE-AGN catalogue could present some contamina¬ 
tion, which does not correlate neither with other surveys 
nor with the GMB, but introduces some bias on the ISW 
amplitude estimation. 

Another interesting aspect is shown in Eig. 8, where 
the correlation between ISW amplitudes is given, for all 
the possible survey-survey combinations. Using our coher¬ 
ent simulations, we have studied the correlation coefficient 
(p) among the estimated A amplitudes from each of the 
six surveys. On each panel of this figure, we show a scat¬ 
ter plot obtained from 1000 simulations, confronting the 
ISW estimation for two surveys. We also plot (red circle) 
the value corresponding to the data. Einally, the correlation 
coefficient is also given. 

As previously mentioned, the CMB-Kappa cross¬ 
correlation is one of the most robust results, since it repre¬ 
sents a detection of the ISW effect, fully obtained from 
Planck data, which simplifies consideration of possible 
sources of systematics present in galaxy catalogues. This 
correlation is nothing but an estimation of the ISW-lensing 
bispectrum (Lewis et al. 2011) induced by the lensing effect 
suffered from the CMB photons, as they pass through the 
gravitational potential. 

This bispectrum represents a bias when determining 
some primordial bispectrum shapes, and one needs to ac¬ 
count properly for it. Within the Planck collaboration, 
this alternative way of measuring the ISW-lensing corre¬ 
lation is carried out in Planck Collaboration XVII (2015). 
In fact, a similar cross-correlation to the one described 
here is also performed in the Planck lensing paper (Planck 
Collaboration XV 2015), as a check to establish the reli¬ 
ability of the lensing map reconstruction. In the following 
section, we summarize these alternative ISW-lensing esti¬ 
mations performed within the present 2015 release. 

3.3. Results on the ISW-lensing bispectrum 

The Planck 2013 results (Planck Collaboration XIX 2014; 
Planck Collaboration XVII 2014; Planck Collaboration 
XXIV 2014) showed for the first time evidence of the 


lensing-ISW CMB bispectrum by using the Planck 2013 
temperature-only data release. The lensing-ISW non- 
Gaussian signal is an independent and direct probe of the 
influence of dark energy on the evolution of structure in the 
Universe, which only relies on CMB data. 

The lensing potential 0 and the CMB temperature T 
are correlated, since it is the same gravitational matter 
distribution at redshifts less than about 2 that leads to 
both the gravitational lensing of the CMB and the ISW 
effect. Moreover, since the gravitational lensing leads to 
changes in the small-scale power of the CMB, and the 
ISW effect affects the large-scale CMB temperature, we ob¬ 
tain a non-zero lensing-ISW bispectrum of a predominantly 
squeezed shape, correlating one large scale with two much 
smaller scales (see, e.g., Goldberg & Spergel 1999; Seljak 
& Zaldarriaga 1999; Hu 2000; Hu & Okamoto 2002; Verde 
& Spergel 2002; Giovi et al. 2003; Okamoto & Hu 2003; 
Giovi & Baccigalupi 2005; Lewis & Challinor 2006; Serra & 
Cooray 2008; Mangilli & Verde 2009; Hanson et al. 2009, 
2010; Smith & Zaldarriaga 2011; Lewis et al. 2011). 

The 2015 Planck release offers us the possibility of in¬ 
cluding polarization in the estimation of the ISW-lensing 
bispectrum. As shown in Cooray & Melchiorri (2006), the 
direct ISW-lensing correlation in E polarization due to re¬ 
scattering of the temperature quadrupole generated by the 
ISW effect is negligible. However, as explained in Lewis 
et al. (2011), there is an important correlation between the 
lensing potential and the large-scale E polarization gener¬ 
ated by scattering at reionization. Because the lensing po¬ 
tential is highly correlated with the ISW signal, this does 
in the end also lead to a non-zero ISW-lensing bispectrum 
in polarization. Although the current high-pass filtering of 
the polarization data reduces this cross-correlation some¬ 
what, it is in principle still detectable. Explicit expressions 
for the ISW-lensing bispectrum template can be found in 
Planck Collaboration XVH (2015). 

In this section we summarize the ISW-lensing estima¬ 
tions performed in three different papers of the present 
Planck 2015 release, and we comment on their compar¬ 
ison. Eirst, as explained in the previous subsection, we 
have implemented an estimator (see Eq. 18) of the ISW- 
lensing bispectrum in terms of the CMB and lensing cross¬ 
correlation (Lewis et al. 2011). An independent implemen¬ 
tation of the same estimator can be found in the Planck 
lensing paper (Planck Collaboration XV 2015). However, 
whereas the latter uses the EEP8 simulations (which in¬ 
clude the actual non-Gaussian signal induced by the lensing 
of the CMB anisotropies), the implementation performed in 
this paper uses Gaussian simulations that form part of the 
set of coherent CMB and LSS tracer maps. Despite this 
difference, both implementations yield very similar results, 
the ISW paper estimator gives A = 1.06 ± 0.33, whereas 
the lensing paper estimator finds A = 0.90 ± 0.28 (both for 
SMICA). 

The Planck primordial non-Gaussianity paper (Planck 
Collaboration XVH 2015) studies the ISW-lensing signal 
primarily to determine the bias this induces on the different 
primordial bispectrum shapes. However, it also gives results 
for the actual amplitude of the ISW-lensing signal. Three 
different estimators have been considered in this paper: the 
KSW estimator (Komatsu et al. 2003); the modal estima¬ 
tor (Eergusson et al. 2010), and the binned bispectrum esti¬ 
mator (Bucher et al. 2010). Whereas the two former meth¬ 
ods are only implemented to work with temperature data. 
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Fig. 8. Correlation among the estimated Aisw amplitudes for the different surveys. The small blue dots are the ampli¬ 
tudes estimated from the simulations described in Sect. 2.3, whereas the large red dot stands for the amplitudes estimated 
from the data. For each pair, the correlation coefficient is indicated. 


the binned bispectrum estimator is also able to include po¬ 
larization. All these estimators use the FFP8 simulations 
to characterize the expected signal and the uncertainties. 
The binned bispectrum estimator finds an amplitude of the 
ISW-lensing bispectrum of A = 0.82 ± 0.27 for the SMICA 
map using both temperature and polarization. The values 
obtained using temperature alone with the KSW, modal 
and binned bispectrum estimators are, A = 0.79 ± 0.28, 
A = 0.72 ± 0.26, and A = 0.59 ± 0.33, respectively. 

We performed a study on a set of 100 FFP8 simula¬ 
tions that have passed through the SMICA component sep¬ 
aration pipeline to investigate any biases in the estimators 
and their correlations. For reasons explained above, the two 
implementations of the estimator based on the lensing re¬ 
construction are not exactly the same. There are also differ¬ 
ences between the three bispectrum estimators; the KSW 
estimator implements the ISW-lensing template exactly 
(since it is separable), while the modal and binned estima¬ 
tors use approximations. Unlike all other templates stud¬ 
ied in Planck Collaboration XVII (2015), the ISW-lensing 
template is difficult to bin and the correlation between the 


exact and binned template is relatively low. Another dif¬ 
ference is that the bispectrum estimators use i^nin = 40 
in polarization, while the lensing reconstruction estimators 
use .^min = Foi* ^11 thcsc reasons we do not expect the 
correlation between the different estimators to be perfect, 
which leads to slight differences in the results. The result 
of the study is that all the bispectrum estimators agree on 
the average value, which is slightly low at around 0.85. The 
other estimators find higher values. Since these same sim¬ 
ulations are used to determine uncertainties on the final 
result, all the error bars have been divided by the average 
that each estimator finds. Regarding the correlations, we 
find that the KSW and modal estimator are correlated at 
about 95 %, while their correlation with the binned estima¬ 
tor is about 80 %. The two lensing reconstruction estimators 
are also correlated at about 80 %, while the correlation be¬ 
tween the two types of estimator classes is about 60-70%. 

Despite these differences, we conclude that all results 
are consistent with the expected value for the ISW-lensing 
bispectrum amplitude A = 1, and that the absence of any 
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ISW-lensing signal {A = 0) is excluded at the level of about 
3(7. 

3.4. Derived cosmological constraints on dark energy 

We have explored the possibility of constraining some cos¬ 
mological parameters through the ISW detection reported 
in 3.2. In principle, the ISW effect depends on the full pa¬ 
rameter set of a dark energy (or curvature) cosmology, but 
the weak overall significance of the signal makes it neces¬ 
sary to restrict parameter measurements to a single or at 
most a pair of parameters, while the remaining parame¬ 
ters need to be constrained from other observations. A ex¬ 
haustive study on dark energy constraints can be found in 
the Planck dark energy and modify gravity paper (Planck 
Collaboration XIV 2015). 

We have assumed a Gaussian shape for the ISW likeli¬ 
hood >C(0), where 0 stands for a general set of cosmological 
parameters: 

-21n[£(e)]=x"(e)-xLn- (19) 

Here the corresponding quadratic x^(0)-functional is given 
by 



Fig. 9. Conditional probabilities on (left panel) and w 
(right panel) derived from the ISW likelihood, based on the 
CMB-NVSS and the CMB-Kappa cross-correlations. 

at 68% CL. These conditional probabilities are shown in 
Fig. 9. 

4. Stacking of CMB temperature and polarization 
data 


x"(e) = 


c, 


TG,obs 




■y-1 


C^G,obs _ 


( 20 ) 

The covariance matrix describes the Gaussian variation 
of the measured spectrum around the theoretical 

expectation Cj^ ^ and is estimated from the 10 000 simula¬ 
tions described in Sec. 2.3 and used in the cross-correlation 
analyses in the previous section: 


= ( 21 ) 

i=l 

where A(Q^) = Cj^{Qo) — Cj^ are the esti¬ 

mates for every single realization i, and is their theo¬ 
retical value. Non-zero off-diagonal entries describe correla¬ 
tions between different multipoles due to broken homogene¬ 
ity, which is mainly caused by masking of emission from the 
Milky Way. 

We expect that the covariance matrix does not change 
strongly with the cosmological model and therefore, that 
the fiducial model 0o given in Sec. 2.3 provides a suitable 
uncertainty characterization for all considered cases. The 
Gaussian likelihood adopted above is the common choice for 
this effect (see e.g., Nolta et al. 2004; Vielva et al. 2006; Ho 
et al. 2008). In our case, this likelihood is used to explore the 
conditional probability of a given cosmological parameter 
(e.g., Oa), keeping constant the remaining cosmology. In 
this case, it is trivial to prove that the estimator is unbiased. 

For simplicity, the data used for the ISW likelihood is 
the joint cross-correlation of the Planck GMB map with the 
NVSS and the Kappa tracers, which already captures 95 % 
of the total detection of the ISW effect (see Table 2). First, 
we have determined the conditional probability (where the 
rest of the cosmological parameters are fixed to the Planck 
fiducial model, Planck Collaboration XIII 2015) for (2 a, ob¬ 
taining the best-fit for (2 a = 0.67 and 0.49 < (2 a < 0.78 
at 68% CL. In particular, we find (2 a > 0 at more than 
3 (7. Second, we have estimated the conditional probability 
on the equation of state parameter of the dark energy, ob¬ 
taining the best-fit for w = —1.01 and —4.45 < re < —1.07 


As an alternative approach to the detection of the ISW 
signal, we can focus on the objects expected to yield the 
strongest effect, namely the largest (tens to hundreds of 
Mpc) voids and clusters in the Universe. In order to mea¬ 
sure the effect produced by individual structures, one can 
stack patches of the CMB anisotropy map centred at the 
locations of superstructures on the sky. Such a stacking 
technique allows us to detect and characterize a signal that, 
otherwise, would be undetectable due to the weakness of the 
ISW effect compared to the primordial CMB anisotropies. 

Following this approach, (Granett et al. 2008a, here¬ 
after GR08) found a potentially significant ISW signal by 
studying 100 superstructures identified in the SDSS DR6 
LRG catalogue. The presence of this signal has since been 
confirmed and more precisely studied with the latest CMB 
data (Planck Collaboration XIX 2014). However, the statis¬ 
tical significance of this detection is still debated, as well as 
its supposed ISW nature (Hernandez-Monteagudo & Smith 
2013; Ilic et al. 2013; Planck Collaboration XIX 2014) and 
the compatibility of its high amplitude with ACDM predic¬ 
tions of the ISW effect from such structures (Granett et al. 
2008a; Hernandez-Monteagudo & Smith 2013; Cai et al. 
2014; Hotchkiss et al. 2014). Moreover, more recent cata¬ 
logues of superstructures have since been used for similar 
studies (Planck Collaboration XIX 2014; Kovacs & Granett 
2015), but none of them has yielded a signal with the same 
level of significance as the GR08 catalogue. 

A crucial point in stacking studies is to determine what 
fraction of the signal detected using this method is either 
due to the ISW effect of the observed structures, or ran¬ 
dom and fortuitous anisotropies of the primordial CMB, 
or a mixture or both. In the present section, we attempt 
to address this question for the results obtained with the 
GR08 catalogue, sine it is to date the only result to appar¬ 
ently show a significant discrepancy with respect to ACDM 
expectations. The main novelty of the present analysis com¬ 
pared to previous works in the literature will be the use of 
a variety of statistical tests that rely on the latest polariza¬ 
tion data from Planck. Indeed, the CMB polarization map 
should prove to be a valuable asset for our purposes; any 
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ISW signal found in temperature is expected to have no 
counterpart in CMB polarization, whereas we expect that 
a primordial CMB signal will be correlated at some level 
with the CMB polarization. Therefore, and despite the lack 
of the largest scales (see Sect. 2.1) in the polarization data, 
it can be used as a discriminant to separate genuine ISW 
detections from false positives due to random primordial 
anisotropies. 

In practice, our objective here will be to answer the 
following questions. Can polarization data help us to prove 
the ISW nature of the GR08 signal? Or disprove it - i.e., 
show that it is actually caused (partially or entirely) by 
the primordial part of the CMB? We should keep in mind 
that the answers to these two questions could very well be 
negative, if the discriminating power of polarization data 
proves to be insufficient for stacking studies. In addition, 
the validity of the GR08 catalogue as an LSS tracer is also 
addressed by stacking patches from the Planck lensing map. 

4.1. Stacking methodology in polarization 

The main procedure for stacking of CMB patches in the 
ISW context has been detailed in Planck Collaboration XIX 
(2014). However, the process for stacking patches of polar¬ 
ization data is not as straightforward as for scalar signals 
like the CMB temperature or the E mode polarization; in¬ 
deed, the Q and U tensorial components are referred to a 
local frame, and patches at different locations cannot be 
directly stacked together. Here instead, we employ a con¬ 
figuration of the Stokes parameters that allows for superpo¬ 
sition; more precisely, we use the following locally defined 
rotation of the Stokes parameters: 

Qr(0) = -Q (0) COS (20) - (0) sin (2(/)) , . 

U,{d) = Q (6>) sin (20) - [/(6>) cos (20) ’ 

where 6 = 9 (cos 0, sin 0) and 0 is the angle defined by 
the line that connects the location considered at the centre 
of the reference system and a position at an angular dis¬ 
tance 0 from the centre. This definition, first proposed by 
Kamionkowski & Loeb (1997), decomposes the linear po¬ 
larization into a radial (Qr > 0) and a tangential (Qr < 0) 
contribution around the reference positions. Komatsu et al. 
(2011) provided a recipe to compute the theoretical T, Qr 
and Uy angular profiles from stacked patches centred on 
temperature peaks, making explicit its dependence on the 
correlations of the CMB primordial anisotropies. 

In practice, we also remove the monopole and dipole 
from the temperature maps outside the mask, before 
computing the signal around each location of the 
GR08 structures. Similarly to the work done in Planck 
Collaboration XIX (2014), we then compute two types of 
profiles from each T, Qr and patch. On the one hand, 
the radial angular profile is obtained as the mean of the pix¬ 
els in rings of fixed width. We choose 150 different angular 
scales from 0° and 15°, with a width of 60 = 0?5 for each 
ring. On the other hand, the value of the photometry profile 
at a given scale is defined as the difference between the av¬ 
erage signal within the disk of radius 0 and the surrounding 
ring of equal area (i.e., between radius 0 and \/26>). In this 
case, 150 angular scales are also taken into account, defin¬ 
ing aperture sizes between 0° and 15°/ a/2 10?6. The 

final step is to compute the average of all profiles (radial or 
photometric) for all of the selected locations. 
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Fig. 10. Stacked patches of the scalar components from 
the SEVEM solution, at the supercluster (first column) and 
supervoid (second column) positions from GR08. From top 
to bottom: T, Te-c^ Te-u^^Ol components. Temperature 
maps are given in pK units. 


Complementary to the temperature analysis, the 
stacked profiles are also computed for the F^correlated 
{Te-c) and the F^uncorrelated (Te-u) temperature maps 
(see Sect. 2.1), as well as for the Planck lensing map, where 
we search for a counterpart to the anomalous temperature 
signal. The stacked images for every map and set of struc- 
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Clusters Voids Total 



Fig. 12. Photometry profiles of the stacked temperature 
patches at the supercluster (first panel) or supervoid (sec¬ 
ond panel) positions from the GR08 catalogue. The third 
panel shows the difference between cluster and void profiles. 
Coloured lines correspond to the different component sep¬ 
aration methods: COMMANDER (red); NILC (orange); SEVEM 
(green); and SMICA (blue). Notice that the four lines are 
almost exactly on top of each other. The dotted black lines 
correspond to the mean values of the null profiles, i.e., com¬ 
puted at the same locations as the real superstructures, 
but in 1000 FFP8 simulations processed through the SEVEM 
pipeline. The shaded regions show the Tier and ±2a uncer¬ 
tainties. Similar levels are obtained for the different com¬ 
ponent separation methods. 


tures considered in this section are shown in Figs. 10 and 

11 . 

4.2. Temperature analysis 

In order to confirm the result presented in Planck 
Collaboration XIX (2014), we carry out the stacking of tem¬ 
perature patches at the locations of the CR08 structures. 


For the whole analysis, we use HEALPix maps at Ngide = 512 
with a filter of FWHM = 20 aremin. The mean radial 
and photometry profiles are computed from each set of 50 
superclusters and supervoids, respectively. Simultaneously, 
we perform the same analysis on 1000 FFP8 simulations 
of CMB temperature, and derive the statistical properties 
of the resulting profiles (mean and standard deviation at 
all scales). Finally, we determine if the profiles measured 
on real CMB data present any significant deviation from 
those derived from the simulations. Since we expect the 
simulated maps to have no correlation with the actual large- 
scale structures of the Universe, this procedure corresponds 
to carrying out a null hypothesis test. 

As we show in Fig. 12, we observe the peculiar shape 
for the profiles already detected in Planck Collaboration 
XIX (2014) using the CMB temperature maps supplied by 
the different component separation methods, i.e., an excess 
of temperature signal at scales around 5° in the photome¬ 
try profiles computed on the supercluster positions, and a 
deficit at scales around 4° in the corresponding supervoid 
locations. The deviation is even more evident if the total 
photometry profiles are computed as the difference between 
the profiles from clusters and voids, as is shown in the third 
panel of Fig. 12. 

A multi-frequency analysis on SEVEM maps is also per¬ 
formed to check if these deviations are monochromatic or, 
conversely, show a specific frequency dependence. In Fig. 
13, we show the mean temperature profiles computed in the 
100 CHz, 143 GHz, and 217 GHz maps. The error bars are 
estimated as the dispersion of the mean profiles computed 
at the GR08 positions in 1000 FFP8 simulations processed 
through the corresponding SEVEM pipeline. As we show in 
Fig. 13, the temperature signal is frequency independent, 
as already checked in Planck Collaboration XIX (2014). 
The uncertainties plotted in the panel correspond to the 
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Fig. 11. Stacked patches of the Qr (top) and Ur (bottom) 
components from the SEVEM solution in /iK, at the super¬ 
cluster (first column) and supervoid (second column) posi¬ 
tions from GR08. 
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Fig. 13. Mean photometry profiles of the stacked temper¬ 
ature patches at the supercluster (first panel) and super¬ 
void (second panel) positions of the GR08 catalogue. The 
CMB data used are the 100 GHz (red), 143 GHz (blue), 
and 217 GHz (green) cleaned maps supplied by SEVEM. The 
dotted black line and shaded regions show the mean, Tier 
and ±2(7 uncertainties of the null profiles, i.e., computed at 
the same locations as the real superstructures but in 1000 
FFP8 simulations processed through the 143 GHz SEVEM 
pipeline. 


143 GHz case only, since the level of the corresponding un¬ 
certainties for the other frequencies is similar. 

In summary, the analyses performed in this section con¬ 
firm that, as expected, the Planck 2015 temperature data 
also exhibit an anomalous signal that can be associated 
with the GR08 catalogue. 

4.3. Polarization analysis 

One of the most attractive developments for this Planck 
release is the possibility of exploring the counterpart in po¬ 
larization of these temperature anomalies. Our motivation 
here is the following: if the 4 a signal measured in tem¬ 
perature is dominated by the primordial part of the GMB, 
it is reasonable to assume that it will have some form of 
counterpart in polarization (although it might not be de¬ 
tectable within the total polarization signal). Gonversely, 
if these temperature deviations are created as a result of 
the presence of large clusters and voids, then no correlated 
signal is expected in polarization. 

In the following, we take several approaches in order to 
make the most of the potentially discriminating power of 
the polarization data, focusing on the study of the GR08 
results. It should be noted, however, that the strength of 
these tests could be diminished by the high-pass filtering of 
the Planck 2015 polarization data release. 

4.3.1. Qr/Ur profile significance estimation 

Since the Qr signal is proportional to the correlation be¬ 
tween temperature and F^mode polarization (e.g., Komatsu 
et al. 2011), it represents a valuable observable for studying 
a potential polarization counterpart to the previously ob¬ 
served temperature signal. It should be noted that, a priori, 
no signal is expected in the 17^ map since it would depend 
on TH correlations, which are null in the standard model. 


The aperture photometry profiles are shown in Fig. 14 
and overall do not present any significant signal at large an¬ 
gular scales (greater than 1°). A notable exception comes 
from the 17^ photometry profile for the voids, which does 
show two significant excesses around 3° and 6°. However, 
since no TH correlations are expected either for the pri¬ 
mordial GMB or for the ISW effect, these features are most 
likely caused either by a fortuitous signal, and/or systemat- 
ics in the polarization map that remain to be characterized, 
which is not in the scope of this paper. The deviations seen 
at angular scales below 1°, especially in the supercluster 
case, are somewhat reminiscent of the expected primordial 
Qr peaks, which appear due to the dynamics of the pho¬ 
ton flows around over-dense and under-dense regions at the 
last-scattering surface (i.e., hot and cold spots in the pri¬ 
mordial GMB, see, for instance Hu & White 1997; Komatsu 
et al. 2011; Planck Gollaboration XVI 2015, for a descrip¬ 
tion of these dynamics). However, this similarity is most 
likely fortuitous since the shape of the temperature pro¬ 
files does not bear much similarity to that obtained from 
the stacking of the extrema of the primordial GMB. In ad¬ 
dition, the GR08 positions do not appear to correspond 
with the positions of GMB extrema. The total polarization 
photometry profiles, computed as the differences between 
clusters and voids, are shown in the third column of Fig. 14. 

We should note, however, that the significance and in¬ 
terpretation of these results are complicated by the use of 
high-pass filtering in the polarization data, which could 
mitigate the signal, in principle, at all scales in the pro¬ 
files (since they include contributions from a large range of 
multipoles). 

4.3.2. Covariance analysis 

We have also tried a different, more general approach to the 
problem by focusing on the following question: what should 
we expect in the Qr/P7x stacking signal, if the GR08 temper¬ 
ature signal originates purely from primordial anisotropies? 

To answer this question, we perform a set of 100 000 
simulations of GMB T, Q, and U maps, using the Planck 
best-fit cosmological model as input. For each one of these 
sets of maps, we derive the T, Qr, and stacked images 
corresponding to the 50 sky positions of the 50 voids of the 
GR08 catalogue (in order to keep the same, potentially rel¬ 
evant, configuration of positions in the sky). For these im¬ 
ages, we derive the radial and photometry profiles, and end 
up with a collection of 2 (temperature and photometry) x 3 
(T, Qr, Ur) X 100 000 profiles. More precisely, for each one 
of the 100 000 sets of GMB T, Q ,and U maps, we can con¬ 
struct the corresponding vector in which we put end to end 
the three radial profiles (of the T, Qr, Ur stacked images) 
and the three photometry profiles. Using these 100 000 vec¬ 
tors, we construct a covariance matrix which contain 
the covariance between any combination of angular scales of 
any of the T/Qr/Ur profiles. We also derive the correlation 
matrix Nij , defined as Nij = Mij / y^^UiMjj. 

In the resulting matrices, we look for the existence of 
significant correlations between a temperature signal, with 
features similar to the GR08 one (i.e., peaking around a 
scale of 4°), and a polarization signal at any scale. The idea 
here is that in the simulated maps and associated stacked 
images that we use here, we can be certain that any stacked 
signal that arises in temperature is fortuitous and due to 
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primordial anisotropies. Starting from this point of view, 
the covariance analysis allows us to obtain a general pic¬ 
ture of how a primordial stacked signal in temperature is 
correlated to its (potential) polarization counterpart. This 
provides us with valuable insight when trying to test the 
hypothesis that the GR08 signal is purely (or partially) pri¬ 
mordial. 

After performing this analysis, the covari¬ 
ance/correlation matrices obtained show that the 
temperature photometry at around 4° is indeed cor¬ 
related with a polarization signal, both in the Qr radial 
and photometry profiles (with a maximum correlation 
around 4° for both). The existence of these correlations is 
quite robust, thanks to the large number of simulations, 
and confirms that if a significant, primordial CMB signal 
appears in temperature, it will have a counterpart in 
polarization. However, it shows that for a GR08-like signal 
only due to primordial GMB, the biggest correlation 
factors with polarization are below about 15 %. Therefore 
a 3 cr signal in temperature would only translate on average 
into 0 . 5(7 signal in Qr, making it effectively impossible to 
detect among the rest of the polarization signal. 

4.3.3. Te-c and Te-u maps 

An alternative to the use of the Qr and components is 
to perform the stacking at the locations of the GR08 struc¬ 
tures, but using the T^correlated and T^uncorrelated tem¬ 
perature maps of the GMB described earlier (see Sect. 2.1). 
If we were to find that most of the GR08 signal is contained 


in the stacked image associated with the E-correlated tem¬ 
perature map, this would be a strong argument towards a 
primordial nature of this signal. Gonversely, if it is found 
mostly in the Euncorrelated temperature stacked image, it 
would give credence to the ISW signal hypothesis. 

In Fig. 15, we show the mean photometry profiles com¬ 
puted from the two aforementioned maps. It appears quite 
clearly that for both the superclusters and the super¬ 
voids, most of the signal originally observed in temperature 
is contained in the E-uncorrelated map, thus apparently 
strengthening the hypothesis of the ISW nature of the sig¬ 
nal, however we have to bear in mind that, by construction, 
the Te-u iiiap contains more power than the Te-c map. 
We remark that, although the polarization data have been 
high-pass filtered, most of the relevant scales responsible 
for the anomalous temperature signal could be still present 
in the analysed maps, since we have checked that the pho¬ 
tometry profiles on high-pass filtered temperature data are 
very similar to those plotted in Fig. 12. On the other hand, 
the ^correlated part of the signal appears to sit within the 
expected values from simulations. 

4.4. Planck lensing convergence map 

Another physical observable explored in the present anal¬ 
ysis is the Planck lensing convergence map. As explained 
in Section 2.2, the Kappa map is proportional to the grav¬ 
itational field and is therefore expected to be correlated 
with the distribution of the large-scale structures, as well 
as the individual objects that generate the ISW signal in the 



02468 10 2468 10 2468 10 

deg deg deg 


Fig. 14. Mean photometry profiles of the Qr (first row) and 
Ur (second row) components stacked at the 50 superclus¬ 
ter (first column) and 50 supervoid (second column) posi¬ 
tions of GR08. The GMB maps used come from COMMANDER 
(red); NILC (orange); SEVEM (green); and SMICA (blue). The 
third column shows the difference between the supercluster 
and the supervoid profiles. The dotted black lines show the 
mean value of our estimator computed with 1000 FFP8 sim¬ 
ulations processed through the SEVEM pipeline at the same 
locations as the real superstructures. Shaded regions show 
the ±I (7 and ±2 a uncertainties of these null profiles; those 
computed for the rest of component separation methods 
reach a similar level. 


20 










Planck Collaboration: The ISW effect with Planck 


CMB. Although complicate projection (and possibly can¬ 
cellation) effects are expected to be involved here, we can 
expect that the stacking of the leasing map at the locations 
of the GR08 structures will give a significant signal with 
respect to the null hypothesis. On the other hand, an ab¬ 
sence of signal could indicate a problem with the structures 
and it could put into question the method and data used 
to identify them in the SDSS and therefore even question 
their existence, or at the very least their reported proper¬ 
ties (sizes, redshifts, etc). The legitimacy of such questions 
is reinforced by recent studies (see Kovacs & Granett 2015) 
that failed to detect some of the GR08 structures in newer 
SDSS data that cover the same survey volume. 

The photometry profiles computed from the Kappa map 
at the positions of the GR08 structures are shown in Fig. 
16. The error bars are estimated with simulations gener¬ 
ated according to the leasing model. We should be cau¬ 
tious here when drawing any conclusions, since the leasing 
map is known to be very noisy. However, it should be noted 
that the photometry profiles for the clusters and voids show 
some relatively significant features, with opposite signs, as 
expected if these profiles result from the averaging of several 
gravitational wells and hills, respectively. The significance 
of these features reaches as high as 3 a for the void profiles. 


although it is hard to pinpoint any typical scale for either 
of the two cases. However, the fact that the profiles tend 
to be positive in the case of clusters and negative for voids 
could be pointing in favour of the ISW interpretation of the 
temperature signal observed. 

To shed light on this question beyond these qualitative 
statements, a estimator is used to evaluate the compat¬ 
ibility of a binned version of these profiles with the null hy¬ 
pothesis, taking into account the correlation between differ¬ 
ent scales. In the case of superclusters, the p-value reaches 
a value of 50 %, which could represent an evidence that the 
observed positive trend is not significant (and without any 
further evidence in favour of the ISW interpretation). In 
the case of supervoids, the same probability is around 2 %, 
such that the observed feature could be a hint of a signif¬ 
icant signal with respect to the null hypothesis from the 
GR08 structures. 

4.5. Summary of stacking 

A stacking analysis of the GR08 positions reveals a tem¬ 
perature signature, as in Planck Collaboration XIX (2014), 
which according to the literature compares poorly with the 
ISW predictions. The major deviation appears at about 
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Fig. 15. Mean photometry profiles of the images for the 50 
superclusters (first column) and the 50 supervoid (second 
column) positions of GR08, when stacking in the Te-c (first 
row) and Te-u (second row) maps computed from the CMB 
maps supplied by COMMANDER (red); NILC (orange); SEVEM 
(green); and SMICA (blue). The dotted black line represents 
the null hypothesis computed as the mean value of the pho¬ 
tometry profile at 50 random positions in 1000 FFP8 simu¬ 
lations processed through the SEVEM pipeline, according to 
the noise properties of the CMB data at the GR08 super¬ 
structure locations. The shaded regions show the ±1 a and 
±2cr uncertainties of these profiles, computed as the dis¬ 
persion of the mean photometry profiles of the simulations. 
The corresponding error bars for the different component 
separation methods reach a similar level. 
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4.5° for clusters and 3.5° for voids. An analysis of different 
cleaned frequency maps that SEVEM provides shows that 
the photometry profiles are not frequency dependent. This 
is both compatible with a pure CMB component or an ISW 
signal, and effectively rules out, for instance a hypothetical 
foreground contribution. 

The use of polarization as a discriminant is the main 
novelty of this analysis with respect to the previous Planck 
results. However, the large-scale information of Planck 2015 
polarization is suppressed with a high-pass filter, and the 
conclusions derived from these data should be therefore 
taken with some caution. The Qr photometry profile is re¬ 
vealed to be mostly compatible with the expected signal 
from random positions. The absence of a counterpart in 
polarization is expected for a contribution caused by a sec¬ 
ondary anisotropy. On the other hand, a theoretical covari¬ 
ance analysis shows that a primary temperature anisotropy 
does have a counterpart in polarization, but at such a weak 
level that it would be difficult to detect. 

The analysis of the T^correlated and T^uncorrelated tem¬ 
perature maps at the GR08 locations supplies a comple¬ 
mentary, and supposedly cleaner way to access the potential 
polarization counterpart. We found that the largest part 
of the temperature excess appears in the T^uncorrelated 
component, and is comparable to the signal recovered 
from a high-pass filtered version of the total tempera¬ 
ture map. Moreover, the stacking of the T^correlated maps 
seems compatible with the contribution of random posi¬ 
tions. Although we cannot conclude that the excess is not 
present in the primordial contribution, we do assert that 
it is compatible with a contribution caused by a secondary 
anisotropy, and therefore with an ISW signal. To summa¬ 
rize, we have found some hints that seem to point towards 
an ISW interpretation of the stacked signal observed in 
temperature at the position of the GR08 superstructures. 
However, our analysis of the current high-pass filtered GMB 
polarization maps cannot yet completely confirm nor in¬ 
validate an ISW origin. These analyses have to be fur¬ 


ther explored with the next Planck data release, where the 
GMB polarization is expected to be recovered at all angular 
scales. 

As an additional test, we performed a similar stacking 
of the same positions in the Kappa map. It revealed a sig¬ 
nificant negative signal in the photometry profile associated 
to the voids positions. This most likely confirms the pres¬ 
ence of supervoids (more precisely gravitational potentials) 
at these positions, which helps clearing some doubts about 
the detection and use of such structures, even if the overall 
amplitude remains unexplained. 


5. ISW map recovery 

As we did in Planck Gollaboration XIX (2014), we ap¬ 
ply the linear covariance-based (LGB) filter first introduced 
by Barreiro et al. (2008), and recently extended in Manzotti 
& Dodelson (2014) and Bonavera et al. (2016) to deal with 
several LSS tracers jointly. 

An alternative approach to estimate the ISW tempera¬ 
ture map that does not require the input of any GMB data, 
but only LSS data, is discussed in Appendix A, following a 
simplified version of the method proposed by Kitaura et al. 
(2010) and Jasche & Kitaura (2010), who estimated the 3D 
gravitational potential of a galaxy network given in terms 
of a redshift catalogue. 

The LGB filter is able to combine all the information 
encoded in the GMB and LSS data about the ISW effect, 
in order to recover an actual map of this weak signal. In 
particular, the LGB filter (Barreiro et al. 2008, 2013) was 
originally developed to recover the ISW map by combin¬ 
ing GMB intensity data and one LSS tracer. The method 
has now been extended to deal with any number of LSS 
surveys (Manzotti & Dodelson 2014), as well as to include 
polarization information (Bonavera et al. 2016) . 



deg deg 


Fig. 16. Mean photometry profile of the stacked images 
from the Planck convergence lensing map, at the loca¬ 
tions of the 50 supercluster (left panel) and the 50 super¬ 
void (right panel) from the GR08 catalogue. The dotted 
black line represents the null hypothesis computed as the 
mean value of the photometry profiles at the correspond¬ 
ing GR08 positions in 1000 simulations according to the 
lensing model. The shaded regions show the ±la and ±2 a 
uncertainties of these profiles, computed as the dispersion 
of the mean photometry profiles of the simulations. 
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Fig. 17. Map of the recovered ISW anisotropies (left column) and the corresponding estimated error per pixel (right 
column) obtained from the combination of the Planck SEVEN CMB map with: the NVSS survey (top); the Planck lensing 
map (middle); and both tracers jointly (bottom). The units here are in Kelvin. 


5.1. Methodology 

We briefly describe here the formalism of the extended 
method that will be used in this paper. In order to con¬ 
struct the filter for n given surveys, the covariance matrix 
C{t) between ISW and LSS data is assumed to be known. 
We note that, at each multipole, C is a square matrix of 
order t = n + 1. To simplify the notation, the matrix is 
written such that the first n elements (whose harmonic co¬ 
efficients are given by m)) refer to the auto- and 

cross-spectra involving only the n LSS tracers, while the 
n + 1 element contains the auto- and cross-spectra that 
include the ISW effect being the harmonic coef¬ 

ficients of the CMB intensity map). Through a Cholesky 


decomposition of the covariance matrix, we construct the 
matrix L satisfying C{£) = L(-^)L^(-^). The estimated ISW 
map s{£^m) at each harmonic mode is then given by: 


?(£,m) = 


i=l 


Lit Z^(L 




^tt 
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n ^( 23 ) 


Lit I IIC- 

j=i 


where Cp corresponds to the power spectrum of the CMB 
signal, without including the ISW contribution. To sim- 
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plify the notation, we have dropped the dependence of the 
Cholesky matrix L on £. We note also that although, in prin¬ 
ciple, the inversion of L should be performed only on the 
n X n submatrix, it is equivalent to use of the full matrix, 
since it is triangular. 

The realistic case of incomplete sky coverage or the pres¬ 
ence of Poissonian noise in the surveys can be accommo¬ 
dated in the previous equation. In particular, the contribu¬ 
tion of the Poissonian noise is simply added to the auto¬ 
spectrum of the corresponding survey. For those data with 
partial sky coverage, their corresponding cross- and auto¬ 
spectra are replaced in the filter by its masked version, i.e., 
correlations among different multipoles in the power spec¬ 
tra are integrated with the MASTER algorithm. 

5.2. Results 

We have applied the LCB filter to the Planck CMB tem¬ 
perature map and to different combinations of the surveys 
described in Sect. 2.2. Before applying the filter, the differ¬ 
ent data sets have been masked using an apodized version 
of the masks shown in Sect. 2. The apodization of the masks 
is performed to reduce the spurious correlations introduced 
in the harmonic domain due to incomplete sky coverage. 
To construct the covariance matrix, we have made use of 
the models described in Sect. 2.3. The different auto- and 
cross-spectra are then transformed to their masked versions 
with couplings computed by the MASTER algorithm (Hivon 
et al. 2002). As in previous sections, when using the Kappa 
and WISE maps, a cut for the lowest multipoles is imposed, 
meaning that these surveys do not contribute to the recov¬ 
ered ISW signal for < 8 (for Kappa) and for < 9 (for 
WISE). Eor the CMB intensity map, we consider the Planck 
SEVEM cleaned CMB map, but similar results are expected 
for the other component separation methods. 

To study the contribution of each data set to the final 
ISW map, we have applied the LCB filter to a total of seven 
different combinations of maps. Eor each of these combina¬ 
tions, we consider two different types of masks to study the 
quality of the recovered map, namely the intersection and 
the union mask. The intersection mask only excludes those 
pixels that are masked by all the data sets considered, since 
the method will reconstruct the ISW signal providing there 
is at least one data map available for a given position in 
the sky (although, as one would expect the reconstruction 
error would depend on the number of observations available 
at each pixel). Conversely, the union mask only keeps those 
pixels which are allowed by all the individual masks and, 
therefore, the reconstruction error will be more uniform in 
the region of the sky considered, since the same information 
is available for all pixels. We note that these masks repre¬ 
sent two extreme cases and they are only used to study 
the quality of the reconstruction. To obtain the recovered 
CMB map, all data sets are used after applying their own 
individual masks. The different combinations used to re¬ 
cover the ISW map are given in Table 3, together with the 
sky fraction allowed by the corresponding intersection and 
union masks. 

Eig. 17 shows the reconstruction attained by combining 
the CMB with NVSS, with Kappa and with both surveys 
simultaneously as well, as their corresponding uncertain¬ 
ties per pixel. The intersection mask has been applied in 
each case. The errors are obtained as the average disper¬ 
sion of the input minus the reconstructed ISW obtained 


Table 3. Mean correlation between the input and recon¬ 
structed ISW maps for different combinations of data sets. 



Intersection 

Union 


mask 

mask 


fsky 

P 

fsky 

P 

CMB and NVSS 

0.84 

0.56 

0.56 

0.58 

CMB and Kappa 

0.73 

0.50 

0.61 

0.51 

CMB, NVSS and Kappa 

0.85 

0.58 

0.51 

0.61 

CMB and WISE 

0.70 

0.48 

0.42 

0.49 

CMB and SDSS 

0.69 

0.53 

0.19 

0.60 

CMB and all surveys 

0.85 

0.60 

0.16 

0.67 

All surveys 

0.83 

0.49 

0.17 

0.61 


from 10 000 coherent simulations of the different data sets. 
The first two cases (CMB plus NVSS and CMB plus Kappa) 
were already presented in the Planck 2013 paper, finding 
very similar results to the ones presented here. As one would 
expect, using both tracers jointly with the CMB (bottom 
row) improves the results with respect to the cases where 
only one tracer is used, although the improvement obtained 
by adding the Kappa map is only moderate. This is due, at 
least in part, to the low multipole cut imposed in this tracer, 
which implies that the lowest multipoles are recovered us¬ 
ing only the CMB. The quality of the ISW reconstruction 
can be further quantified by calculating the correlation p 
between the input s and reconstructed s maps using sim¬ 
ulations. Before calculating the correlation, the monopole 
and dipole are subtracted from the input and reconstructed 
maps outside the considered mask. Table 3 gives the aver¬ 
age correlation obtained over 10 000 simulations (outside 
the union and intersection masks) estimated for each sim¬ 
ulation as: 

^ ~ l^s) ~ l^s) ^ (24) 

1/^1 

“ eW 

where the sum runs over all the pixels allowed by the con¬ 
sidered mask and the weights at each pixel uji have been 
estimated from the error map ai shown in the right column 
of Eig. 17. The quantities ag and ag are the dispersion of 
the input and reconstructed map for each simulation ob¬ 
tained with the same weights, while ps and pg correspond 
to the weighted mean values of the same maps. Eor the 
union mask, an average correlation coefficient of 0.61 is 
found when NVSS, Kappa and the CMB are combined, 
to be compared to the cases when only one tracer is used, 
i.e., 0.58 and 0.51 for NVSS and Kappa, respectively. 

Eig. 18 gives the ISW signal reconstructed from CMB 
and the WISE surveys (top) and from CMB and the SDSS 
surveys (bottom) as well as their corresponding uncertain¬ 
ties. A bright red area is seen in the northern Galactic re¬ 
gion, just above the central part of the mask, which can be 
identified with systematics present in the WISE catalogues 
(see Eig. 4). Due to the cut at low multipoles imposed in 
these surveys, the structure at the largest scales is sup¬ 
pressed in the reconstruction, which is reflected in a larger 
uncertainty. The correlation between input and reconstruc- 
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Fig. 18. Map of the recovered ISW anisotropies (left column) and the corresponding estimated uncertainty per pixel 
(right column) from the combination of the Planck SEVEM CMB map with: the two WISE surveys (top); and the two 
SDSS tracers (bottom). The units here are Kelvin. 


tion outside the union mask is 0.49, the lowest value found 
among all the considered cases (see Table 3). Regarding the 
reconstruction using CMB and the SDSS surveys (bottom), 
we find a large signal in the relatively small regions observed 
by these surveys. This is refelcted in a mean correlation be¬ 
tween input and reconstruction of 0.60 in the region allowed 
by the common mask, showing that the SDSS provides a 
sensitive tracer of the ISW effect. 

The top row of Fig. 19 shows the reconstructed ISW sig¬ 
nal, as well as the estimated uncertainty, obtained from the 
CMB map together with the six mentioned surveys (NVSS, 
Kappa, the two WISE and the two SDSS surveys). As one 
would expect, by combining all the available information, 
we obtain the best ISW map with a reconstruction uncer¬ 
tainty of around 14 jaK and a mean correlation coefficient 
of 0.67 outside the union mask. This corresponds to a max¬ 
imum S/N ratio greater than in certain regions of the sky. 
If the intersection mask is considered, the correlation co¬ 
efficient is 0.60 obtained over 85 % of the sky. Finally, the 
last case (middle row) gives the ISW map reconstructed us¬ 
ing only the six surveys, without including the CMB, which 
corresponds to the first term of Eq. (23). It is apparent that 
removing the CMB degrades the reconstruction, especially 
in those areas where less surveys are available, and this de¬ 
creases the correlation coefficient to 0.60 (union mask). To 
show the contribution given by the CMB to the recovery of 
the ISW map, the difference between these two reconstruc¬ 
tions is also shown in the bottom panel of the figure. The 


intensity CMB mask has been applied and the monopole 
and dipole removed outside this mask. We note that, as ex¬ 
pected, the structure of this map mainly reflects that of the 
large scales of the intensity CMB data presented in Fig. 1. 

It is interesting to point out that some common struc¬ 
tures are visible among the different reconstructions, al¬ 
though the maps are not expected to look exactly the same, 
since each survey traces the ISW effect in a different way 
and, thus, each considered LSS tracer provides a partial 
reconstruction of the ISW signal. 

As already mentioned, the structure of the error maps 
given in the right columns of Figs. 17, 18, and 19 reflects 
the different sky coverages of the surveys, showing the con¬ 
tribution of each data set to the final ISW reconstruction. 
This can be further explored by comparing these different 
structures with Fig. 20, which shows the intersection re¬ 
gions defined by the CMB and survey masks. Each colour 
corresponds to a region where the intersection of a differ¬ 
ent sets of masks occurs; the dark blue area is observed by 
all data sets whereas the dark red region gives those pixels 
that are not observed by any of the data sets. 

6. Conclusions 

We have presented a study of the ISW effect using the 
Planck 2015 data release, which provides higher sensitiv¬ 
ity temperature anisotropy maps with respect to the previ¬ 
ous Planck 2013 release, as well as CMB polarization data 
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Fig. 19. Maps of the recovered ISW anisotropies (left column) and the corresponding estimated uncertainty per pixel 
(right column) from the combination of the Planck SEVEM CMB map and all the surveys (NVSS, WISE-AGN, WISE- 
GAL, SDSS-CMASS/LOWZ, SDSS-MphG, and Kappa; top), and only considering the information from these LSS tracer 
surveys (middle). The bottom panel gives the difference between both reconstructions, with the GMB intensity mask 
applied. The units here are Kelvin. 


at angular scales below 5°. Gompared to our past publica¬ 
tion (Planck Collaboration XIX 2014), we have extended 
the analysis in the following ways. 

Eirst, we have included additional galaxy (WISE-GAL) 
and AGN (WISE-AGN) catalogues from the WISE sur¬ 
vey as LSS tracers to be correlated with the four Planck 
CMB maps (COMMANDER, NILC, SEVEM, and SMICA). These 
tracers, in combination with the NVSS radio catalogue, 
the photometric luminous galaxy (LG) catalogue from the 
Baryonic Oscillation Spectroscopic Survey (BOSS) of the 
SDSS III (SDSS-CMASS/LOWZ), and the photometrically 


selected galaxies from the SDSS-DR8 catalogue (SDSS- 
MphG), yield a detection of the ISW signal at 2.9 cr. This 
detection is dominated by the NVSS catalogue (2.6a), 
while the combination of the two SDSS catalogues provides 
a 2.7 a level, and the two WISE render a 1.9 a signal-to- 
noise ratio. 

Second, we have also improved the characterization of 
the ISW effect through the ISW-lensing bispectrum, since 
the higher signal-to-noise ratio of the Planck 2015 tem¬ 
perature data and the new polarization data allows us to 
improve the reconstruction of the Planck lensing signal. In 
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particular, we have increased the detection achieved in the 
previous release by about 20%, reaching rougly a 3 a detec¬ 
tion. We have performed a new analysis in which the Planck 
ISW-lensing is combined with the cross-correlation of the 
Planck CMB with all the previously mentioned LSS trac¬ 
ers, obtaining a total detection of the ISW effect at the 4cr 
level. The four CMB maps provide similar detection levels 
for all the cross-correlation combinations. 

Third, we have investigated the anomalous nature of 
the ISW signal detected through the stacking of the 
CMB anisotropies at the positions of known superstruc¬ 
tures (Granett et al. 2008a). We have confirmed that the 
aperture photometry profiles around the 50 supervoids and 
50 superclusters of the GR08 catalogue exhibit a maximum 
amplitude of —11/iK (at scales of around 3.5°) and 9/iK 
(at scales of around 4.5°), respectively. These amplitudes 
are much larger than expected in the context of the stan¬ 
dard AGDM scenario. We have used the Planck polariza¬ 
tion data to further explore the origin of this signal. We 
do not find evidence for a positive correlation of this signal 
in the polarization data, indicating that the origin of the 
temperature signal is, indeed, compatible with a secondary 
anisotropy, as expected for the ISW effect. These aperture 
photometry results are extremely consistent for the four 
CMB polarization maps, as well as for the SEVEM cleaned 
frequency maps at 100, 143, and 217 GHz. This greatly re¬ 
duces the possibility that this signal is significantly affected 
by contamination from residual Galactic and extragalac- 
tic foregrounds. Similar conclusions are obtained through 
the analysis of the T^-correlated and T^-uncorrelated coun¬ 
terparts of the temperature signal; excess is found only in 
the latter, as expected for the ISW effect. Finally, we have 
also stacked patches of the Planck leasing map at the loca¬ 
tions of these superstructures, and find a positive correla¬ 


tion for both clusters and voids, which offers extra evidence 
in favour of the ISW hypothesis. 

Finally, we have improved the recovery of the ISW fluc¬ 
tuations on the sky by using a generalization of the lin¬ 
ear covariance-based filter. In particular, we have used the 
five galaxy catalogues mentioned above and the Planck 
lensing convergence map to infer a map of the secondary 
anisotropies associated with the ISW effect caused by the 
LSS traced by these surveys. Using simulations, we have 
been able to provide an associated rms map with a mean 
value of 14/rK per pixel of about 1°. Our ISW reconstruc¬ 
tion provides regions where the ISW fluctuations are recov¬ 
ered at more than 2 a. We have also explored an alternative 
approach to estimating a map of ISW anisotropies by at¬ 
tempting a direct inversion of the density field as traced 
by the 2MASS photometric redshift catalogue into its cor¬ 
responding gravitational potential field. The typical rms 
of the ISW effect induced by these nearby structures is, 
as expected, very low 0.6 /iK), and this is well below 
the level of the measured large angular GMB fluctuations. 
Nevertheless, the angular power spectrum of the ISW ef¬ 
fect produced by these structures is accurately recovered 
for i < 20. 

Therefore, the cross-correlation of the Planck CMB 
maps with different tracers of the LSS confirms the detec¬ 
tion of the ISW effect at the expected level for the AGDM 
model. The current detection level could be slightly im¬ 
proved, on the CMB side, by analysing the next Planck 
release, which will include large-scale polarization data. In 
addition, the use of the future full polarization Planck data 
could be very important to probe further the nature of the 
ISW stacked signal, in a more complete manner, since the 
current analysis (limited by the high-pass filetring) provides 
only a partial view of the problem. 



Fig. 20. Regions defined by the different intersections of 
the masks considered for the recovery of the ISW signal, 
using the LCB filter. To produce this figure, for each mask 
we construct a map with a constant value given by 
in the excluded pixels and zero otherwise. These maps are 
then added together, producing the pattern seen in the fig¬ 
ure. In particular, we choose km = {2,4,8,16,32,64} cor¬ 
responding to the masks used for Kappa, NVSS, SDSS- 
CMASS/LOWZ, SDSS-MphG, CMB, and WISE, respec¬ 
tively (the square root function is introduced to allow for a 
better visualization). 
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Appendix A: Construction of an ISW map from 3D 
galaxy surveys 

In this Appendix we present an alternative approach to es¬ 
timate the ISW temperature map that does not require the 
input of any CMB data, but only LSS data. We first outline 
the methodology for a generic survey, and then apply it to 
the 2MPZ survey, which characteristics, in particular, the 
uncertainty associated to the galaxy photometric redshift 
estimates, are small enough to have a proper gravitational 
potential reconstruction. In addition, this survey traces the 
very local universe, whose contribution to the total ISW 
signal has been pointed out by some previous works (e.g., 
Francis & Peacock 2010; Rassat et al. 2013) to be related 
to the large-scale Planck anomalies. 


The approach consists of using redshift information in 
galaxy catalogues to provide a full 3D gravitational poten¬ 
tial reconstruction, which, under the assumption of a given 
cosmological framework, can be trivially converted into an 
ISW map estimate, for linearly evolving structures. We re¬ 
mark that in this case very high redshift precision is not 
required, since the gravitational potential sourcing the ISW 
effect is coming from large scales (at or above 100 h~^ Mpc 
typically), leaving room for redshift uncertainties at the 
level of Az = 0.01-0.03. On such large scales, redshift space 
distortions can be safely ignored. 

A.l. Methodology 

The procedure must invert a galaxy density field into a po¬ 
tential field in a given region of the Universe that is limited 
by the selection function of the survey and the sky mask. 
This is done by applying the Poisson equation in Fourier 
space, and expressing the gravitational potential in terms 
of the density contrast, namely 

-kHk = ^H^n^a-\z)5k. (A.l) 

In this equation, stands for the Fourier transform of the 
gravitational potential, Qm denotes the total matter density 
parameter and Hq corresponds to the Hubble constant. The 
factor a~^{z) = 1 z corresponds to the inverse of the 
cosmological scale factor, and 6k is the (time dependent) 
matter density contrast Fourier transform for the mode k, 
as estimated from the galaxy density. The time dependence 
of the gravitational potentials is thus given by these last 
two quantities. The use of the Poisson equation is justified, 
since we are considering scales that, despite being larger 
than typical density clustering lengths, are well inside the 
horizon. When handling the equation above, the presence of 
an effective volume mask (induced either by the sky mask 
and/or the survey selection function) may introduce biases 
in our gravitational potential estimates. In order to handle 
this, we choose to conduct our particular Poissonian data 
augmentation^ consisting of the following steps 

— We place the galaxies in a regular 3D grid in comov¬ 
ing coordinates. For that purpose we use the central 
value of the redshift assigned to each source. As will 
be shown below, we find that, when accounting for all 
other sources of uncertainty, the level of uncertainty as¬ 
sociated with errors in the photometric redshifts is sub¬ 
dominant and thus can be neglected. 

— In those grid cells excluded by the sky mask, we in¬ 
troduce a number of mock galaxies that comes from a 
Poissonian realization with an average galaxy number 
density equal to the average number density of cells at 
the same distance but not excluded by the sky mask. 

— In all grid cells, we introduce another set of ran¬ 
domly, Poisson distributed mock galaxies, in order to 
make the selection function constant with respect to 
depth/redshift. 

We assume that the selection function of the survey 
and the sky mask can be factorized separately, in such a 
way that the selection function depends exclusively on the 
depth/redshift. In practice, this may not be the case for 
regions with high extinction, but we assume that most of 
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these regions should be discarded by the sky mask of the 
survey. 

This procedure should provide a homogenized galaxy 
density field in the entire 3D grid, which can be inverted 
into a gravitational potential field. By conducting a set of 
simulations for each of the stages of our Poissonian data 
augmentation approach, it is possible to assess the depen¬ 
dence of the resulting gravitational potential field on each 
step. This potential reconstruction method is a very simpli¬ 
fied version of more sophisticated approaches of inversion 
of observed galaxy surveys (e.g., Kitaura et al. 2010; Jasche 
& Kitaura 2010). 

Finally, the gravitational potential time derivative is ob¬ 
tained from the 3D gravitational potential after imposing 
that it follows linear theory predictions of the reference cos¬ 
mological model. In other words, we express the derivative 
of the gravitational potential field with respect to the radial 
comoving distance r] according to 


d^k 

df] 




1 + 


_d\ogD^'] 4 
d\og{l + z) k‘^ 


= -^kH{z)a{z) 


^ dlog Ds 
dlog{lPz) 


(A.2) 


In this equation, H{z) denotes the Hubble parameter, Hq 
its current value, and Ds{z) the redshift dependent linear 
growth factor of the density perturbations. The final ISW 
map is obtained after integrating the d^k/dr] 3D grid along 
the line of sight. 


A.2. Results 

Now we present the results of inverting the 2MPZ survey 
into a gravitational potential field whose time derivative is 
then projected along the line of sight. 

We place the galaxies of this survey in a 3D grid of 
128^ cells of 6h~^ Mpc on a side, centred upon the ob¬ 
server. This means that the maximum redshift considered 
is ^max — 0.13, and that more than 85% of the 2MPZ 
sources are actually placed inside the grid. This choice 
of ^max is motivated as a compromise between sampling 
a large cosmological volume and having a representative 
amount of galaxies tracing the potential wells; increasing 
^max degrades potential reconstructions at large distances 
from the observer and provides little information about the 
ISW signal. 

Provided that the ISW effect is generated out to red- 
shifts z ^ our choice of ^^max should contain a small 
fraction of the total ISW signal generated in our visible 
Universe. Nevertheless the contribution of this relatively 
nearby cosmological volume to the low multipole anisotropy 
power has been claimed not to be completely negligible, and 
it has been argued that it may be of relevance in the context 
of the CMB large-angle anomalies (e.g., Francis & Peacock 
2010; Rassat et al. 2013). 

When conducting the reconstruction, we impose a sky 
mask for all pixels with |6gai| < 10° for which the Galaxy 
heavily impacts the selection function of the survey. The 
result of the inversion of the galaxy density field into the 
gravitational potential field and its time derivative is shown 
in Fig. A.I. The recovered ISW map resembles the large- 
scale structure of the projected density map (Fig. 4). The 
positive structure of the ISW map traces, in the North 
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Fig. A.l. Map of the recovered ISW signal from the 2MPZ 
catalogue. Units here are Kelvin. 



Fig. A.2. Comparison of the ISW-recovered angular power 
spectrum from 2MPZ data (black solid circles) with the¬ 
oretical expectations (blue solid line) and the average of 
ideal (green solid line) and realistic (red solid line) simu¬ 
lations of our density-to-potential inversion algorithm (see 
text for details). 


Galactic hemisphere, the presence of well known superclus¬ 
ters like Ursa Major, Virgo, Centaurus, or Hydra, while in 
the southern hemisphere, at slightly negative Galactic lati¬ 
tudes, the most prominent negative spot corresponds to the 
Local Void. 

We find, however, that the amplitude of the recovered 
ISW map is too small to contribute significantly to the to¬ 
tal CMB map on the largest scales. In Fig. A.2 filled black 
circles display the angular power spectrum of the recov¬ 
ered ISW effect from 2MPZ. The rms of this map is dom¬ 
inated by the quadrupole, whose amplitude is found to be 
^ 2 M^pz,isw _ 44 ^Q 37 ytiK)^, driving the rms map to be 

at the level of only 0.56 fiK. Thus our estimated quadrupole 
amplitude of the ISW map generated by 2MPZ seems to be 
in tension with the estimate of Rassat et al. (2013), since 
those authors quote a theoretical expectation for the ISW 
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quadrupole of (12 ±10 ijlKY ^ i.e., almost two orders of mag¬ 
nitude above our estimate. 

We next compare the amplitude of the angular power 
spectrum of our recovered ISW map with theoretical ex¬ 
pectations. For this purpose, we make use of a modified 
Boltzmann code that provides the ISW angular power spec¬ 
trum for a generic galaxy sample that is probing the large- 
scale structure under the same selection function as the one 
estimated for 2MPZ. We remark that this estimate of the 
ISW amplitude is independent of the bias of the galaxy sam¬ 
ple. Such a prediction is provided by the thick blue solid line 
in Fig. A.2. We can see that the recovered ISW power spec¬ 
trum is significantly higher than this expectation. In order 
to understand this, we run 100 Monte Carlo simulations of 
Gaussian density fields in the same 3D spatial grid used 
for the density-to-potential inversion of real data. These 
Gaussian simulations are obtained from a AGDM matter 
power spectrum corresponding to our fiducial cosmological 
model at z = 0. In this set of MG ideal simulations we only 
impose the radial selection function of the 2MPZ survey 
at the time of conducting the line of sight integral of the 
time derivative of the gravitational potentials, but ignore 
all effects of radial selection function, photometric redshift 
errors, and shot noise, when producing the potential maps. 
The green solid line provides the average angular power 
spectrum obtained from this set of simulated maps. The 
agreement of this computation with the theoretical expec¬ 
tation is very good for multipoles < 20; artefacts related 
to the projection of the finite grid cells on the sky introduce 
spurious power that becomes dominant on smaller scales. 

We also run a second set of MG realistic simulations 
which are based upon the same set of Gaussian simulations 
just described above, but after including the impact of the 
2MPZ radial selection function, photometric redshift errors, 
and shot noise, as was required for real data. The photo¬ 
metric redshift errors were simulated by adding a normal 
deviate of rms, = 0.015 to the ’’correct” redshifts of the 
simulated galaxies. We note that in both sets of simulations 
the initial 3D Gaussian matter density field is identical, but 
in this case Poissonian augmentation was required: (I) in 
the sky-mask excluded areas; and (2) at large redshifts, 
in order to avoid radial galaxy density gradients associ¬ 
ated with the radial selection function. The average angu¬ 
lar power spectrum from this set of simulated ISW maps is 
displayed in Fig. A.2 by the red solid line. We find in this 
case much better agreement with the results from the real 
2MPZ catalogue. From this set of simulations we obtain the 
uncertainties for the recovered 2MPZ ISW power spectrum 
multipoles. 

Despite the uncertainties in the amplitude of the recov¬ 
ered ISW map, this analysis shows that it is highly implau¬ 
sible that, in a standard AGDM scenario, the ISW gen¬ 
erated by the gravitational potentials hosting the 2MPZ 
galaxies can significantly modify the large-scale pattern of 
the GMB; the expected quadrupole of the ISW is about 
three orders of magnitude below the total GMB quadrupole 
amplitude. 

These two sets of simulations should provide a fair de¬ 
scription of the total error budget introduced by our ap¬ 
proach. In Fig. A.3, filled red triangles display the ratio 
of the rms of the recovered angular power spectrum mul¬ 
tipoles over their average value for the realistic set of sim¬ 
ulations. For instance, for the recovered quadrupole this 
plot shows that the rms of the quadrupole amounts to 



Fig. A.3. Relative uncertainty in the recovered ISW an¬ 
gular power spectrum multipoles. Red triangles and green 
circles refer to the realistic and ideal sets of MG simulations, 
respectively. 


o 



Fig. A.4. Correlation coefficient between the ISW recov¬ 
ered maps in the realistic and ideal sets of simulations. 

roughly 90% of its amplitude. We note that the quoted 
uncertainty in the recovered angular power spectrum mul- 
tipoles includes the contribution from cosmic variance. For 
the sake of comparison, the solid blue line depicts this ratio 
for the case of a pure Gaussian field without any coupling 
between different multipoles and the same sky coverage as 
for 2MPZ and the Monte Carlo simulations; in the case this 
ratio obeys the simple form \/2//sky/(2/ ± I). This trend is 
closely followed by the output of the ideal MG simulations 
ignoring the impact of the radial selection function and the 
Poissonian augmentation (green filled circles). For multi¬ 
poles below = 20 errors in the recovered angular power 
spectrum multipoles are close to the Gaussian prediction, 
but on smaller scales errors associated with the line of sight 
integral become dominant. 

Finally, in order to provide an estimate of how the ISW 
maps recovered by our technique actually resemble the real. 
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underlying ISW maps, Fig. A.4 displays the correlation co¬ 
efficient between the ISW recovered map under the ideal 
and realistic set of simulations. This correlation coefficient 
is defined as n = that is, 

the ratio of the angular cross-spectrum of each pair of maps 
over the square root of the product of the auto spectra. The 
correlation coefficient is, on average, about 0.8 in the multi¬ 
pole range £ = 2-10, and above 0.70 in = 2-20. On smaller 
angular scales there is little ISW power and spurious power 
erases any ISW information in the recovered maps. 

Finally, we performed a new realistic set of MC sim¬ 
ulations for which the uncertainty in the radial distance 
to galaxies associated with photometric redshift errors is 
switched off. We found very little difference in the uncer¬ 
tainty in the recovered ISW angular power spectrum mul¬ 
tipoles and in the correlation coefficients, suggesting that 
photo- 2 : errors at the level az = 0.015 are not a dominant 
source of uncertainty for our ISW reconstruction approach. 
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